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Introduction 



The work described in this thesis deals with the structure of the nucleon within the context of 
perturbative Quantum Chromodynamics (QCD). We determine a new generation of unpolarized 
dynamical parton distributions using the latest experimental and theoretical advances in the field 
and paying special attention to the implications of their dependence on the different prescriptions 
utilized in their derivation. Most of the results are to a large extent based on references [1-5]. 

While being confined within hadrons, quarks behave in hard scattering processes as almost 
free particles, i.e., they break the scaling expected in the naive parton model only logarithmically. 
A fundamental property of QCD is the ability to accommodate these rather different aspects of the 
strong interaction in a consistent picture. Qualitatively, the strong force gets stronger as quarks 
move apart; the eventual creation of new quarks and their combination into color neutral hadrons 
{confinement) gets favored before reaching the hadronic size. On the other hand, it gets weaker 
as their distance goes to zero {asymptotic freedom), which permits a reliable perturbative descrip- 
tion of short-distance (large momentum transfer) phenomena. Long-range (soft) interactions are 
related to the structure of hadrons and are not to be described perturbatively. 

This separation of scales or factorization means in practice that hard scattering cross-sections 
involving initial hadrons are described as products of factors describing the collisions sequen- 
tially: the interactions happening "before" the collision, i.e., the structure of the incoming hadrons, 
the collision between partons and, eventually, the formation of the detected hadrons. In the case 
of hadronic final states we will only consider inclusive reactions; for this reason, it suffices for 
our purposes to consider the first two factors. 

The interaction between partons is described by hard scattering partonic cross-sections {coef- 
ficient functions) which contain only short-distance interactions and therefore are (in principle) 
calculable in perturbative QCD (pQCD). The structure of incoming hadrons is described by par- 
ton distribution functions (PDFs). Roughly speaking (at lowest order), they represent the prob- 
ability of finding a particular parton (with certain kinematical properties) within the described 
hadron. It should be clear by now that they contain long-range non-perturbative interactions 
and therefore are not calculable in perturbation theory but rather to be extracted from experi- 
ment. They are, however, independent of the scattering process under consideration {universal- 
ity), which means that, once they are determined from a reduced set of data, they can be used to 
calculate observables; in this sense pQCD is a fully predictive theory. 

Although in principle the partonic description of hadrons is completely general (and even 
extendable to other objects, e.g. photons), the most relevant PDFs are by far those of the nucleon. 
They constitute an ineludible ingredient in any calculation involving initial nucleons, for instance. 
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at the Large Hadron Collider (LHC), which will be the center of attention for particle physicists 
over the next years. Their optimal determination is therefore essential for future achievements 
and we concentrate on them in this work. 

Even after three decades of intensive research and the fact that there are several groups work- 
ing on the field, there are still outstanding issues on the determination of nucleon PDFs, which is 
to some extent model dependent and, in any case, different for different analysis. Which sets to 
use for a particular calculation is somehow a matter of taste but the reader should be aware that 
there is a definite correspondence between partonic cross-sections and PDFs sets, so that one has 
to choose PDFs suitable for combination with the particular partonic cross-section to be used. 
This is sometimes ignored or neglected in the literature, leading to theoretically inconsistent (or 
simply wrong) results. 

First of all, since the PDFs are determined from (global) fits to experimental data, the re- 
sults depend on the kind of data (and the particular sets) which are used and the sensibility of 
the distributions to them. Furthermore, the experimental errors of the data induce an intrinsic 
"error" of the resulting PDFs which is now possible to estimate. The distributions are usually 
extracted as parametrizations of contour conditions for a set of differential equations {evolution), 
and the results are therefore biased by the initial parametrization. We use the so-called dynami- 
cal (or radiative) model, which has proved to be a reliable and highly predictive approach to the 
determination of PDFs. 

There are also intrinsic theoretical uncertainties on the determination of the PDFs which are 
translated into the reliability of the predictions of observables. As mentioned, one works in 
perturbation theory, so that there is always an uncertainty of the size of the last terms neglected 
on the perturbative series (and even a certain freedom in whether or not to include some terms 
and effects). In general a leading order (LO) description is not satisfactory, the standard at the 
moment is the next-to-LO (NLO) approximation which works quite well for most observables. 
The (massless) evolution of the PDFs and some important processes can now be calculated up to 
NNLO and some (few) calculations have been carried out to even higher orders. Other sources of 
dependence/uncertainty are the factorization and renormalization prescriptions and the treatment 
of heavy quark masses. We address all these issues through the thesis. 

After introducing basic elements of pQCD and establishing our theoretical framework, we 
discuss to some extent the dynamical model and compare it extensively with the approach to 
PDFs used by most other groups {"standard") by examining the results of analyses carried out 
under identical conditions in both approaches. Furthermore, parton distributions sets (with uncer- 
tainties) at LO, NLO and NNLO of QCD, using different factorization schemes (MS, DIS) and 
different treatments of heavy quark masses (FFNS, VFNS) are extracted and compared. 

In the last part of the thesis we turn to the analysis of some of the applications of our results. 
We briefly outline the astrophysical impUcations of the dynamical predictions before focusing 
on collider phenomenology. There we study to some extent the relevance and perturbative sta- 
bility of the longitudinal structure function of the nucleon and analyze extensively the role of 
heavy quark flavors in high-energy colliders by comparing the predictions obtained for several 
processes utilizing different prescriptions for the treatment of heavy quark masses. In addition, 
we show how isospin violations in the nucleon help to explain the so-called "NuTeV anomaly". 
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1.1 Renormalization and the Running Coupling 

In renormalizable quantum field theories the logarithmic ultraviolet^ divergences stemming from 
virtual-particle loops can be summed to all orders using the renormalization group. This proce- 
dure replaces the (unphysical) hare coupling constant appearing in the formulation of the theory 

with an effective or renormalized coupling constant = 'f^^ which depends on the substraction 
point, an arbitrary finite mass scale known as the renormalization scale iir. 

The independence of physical quantities on iir is reflected in the so-called renormalization 
group equations, which are solved by implicitly introducing a function as{li^) known as the 
running coupling and given by: 

The observables appearing in scattering processes are usually expressed as functions of ratios of 
invariants {scaling variables) and a single scale (/x) with dimensions of energy (typically related 
to the center-of-mass energy s, e.g. }X ^/s). After renormalization all the dependence of ob- 
servables on the energy scale appears through the running coupling in a well prescribed form (for 
more details see, for instance, [6-10]). 

In contrast to other theories, the ^-function is negative for QCD and therefore the coupling 
decreases as the scale increases (asymptotic freedom), which makes a perturbative treatment of 
the strong interactions at high scales plausible. We can then expand the ^-function itself in 
powers series; the miming coupling is obtained from: 

'^=-£/i.«rv-) (1.2) 

where the expansion coefficients depend in general (the first two are independent) on the 
renormalization scheme used. In this work we will use the MS scheme, in which in addition to 
the divergences only ln4;r — Ye is absorbed into the running coupling; the first coefficients are 
given [11-15] by: 

o 11 2„ o _ino 38„ oMS _ 2857 5033„ . 325„2 /i ox 

P0 = 11— 3«R, pi = lUZ— yWR, P2 — ^ 18~"«+54:"R C^-j) 

being the appropriate number of active quark flavors at the scale n, which is determined by 
convention. 



'Related to very high (— » oo) energy/momentum or equivalently short (— » 0) distance phenomena. 
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In principle, it is possible to fix hr to the number of light flavors and treat the heavy quark 
contributions as separated corrections of the form In ^ , being m the mass of a generic heavy 
quark. However, the most common approach in the MS scheme is the use of a variable number of 
active flavors taken to be the number of flavors which mass satisfies m^<n^. The jS^ coefficients 
are then discontinuous at the thresholds ix=m and the solutions for at both sides of the thresh- 
olds are subject to adequate matching conditions [16-18]. Up to two loops (m= 1) the ruiming 
coupling is continuous; at three loops a marginal term ( — ^a,) appears, although its effects 
are numerically insignificant (some few per mil difference in the relevant regions) and we have 
neglected it in our NNLO analyses. Using this variable number of active flavors the contributions 
arising from heavy quarks are automatically included (resummed) and consequently the stability 
of the perturbative expansion is improved. 

Turning now to the integration of Eq. 1 .2, it is conventional to introduce a parameter Aqcd 
(henceforth just A) to be extracted from data, i.e., determined by a contour condition at a reference 
value, often chosen to be as{M^). A depends in general on the order considered in the expansion, 
the renormalization scheme used, the number of active flavors and the explicit form of the solution 
adopted. Up to next-to-leading-order (NLO) a straightforward solution is: 



ln^ = — -^ln(— + ^) (1.4) 



where the first term alone is a solution of the leading order (LO) equation^. This equation consti- 
tutes the definition of the A parameter since it is determined from Eq. 1 .2 up to a multiplicative 
constant. In the same way, a solution of the NNLO (m = 2) equation may be found by direct 
integration: 

to^^ ' -4l„f ' + ft +fcU ^ftfe-ft^ a,c.a„/ ^ft^ + ft 
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(1.5) 

Note that for = 6 the argument of the square root is negative and this formula does not give a^; 
in this case a solution is given by: 

+i32-j3oi32-i3ii32A 




Since Eqs. 1.4 to 1.6 constitute implicit solutions for a^, a numerical iteration is required for 
the exact determination of the coupling beyond LO and one could equally integrate numerically 
Eq. 1.2 directly, without introducing the A parameter. In principle one could also expand as{li^) 
in a power series in as{Hg) for some fixed fig-, up to ^(a^) we get: 

fl,(M') = asill^g) - a^.illi) ft In g + a^illi) (^P^ln^ g - ft In |) , (1.7) 



^of course A'^ / A'^", Oj =aj(/j^), etc.; we try to keep the notation as compact as possible 
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however, the uncertainty introduced by higher order terms makes these kind of expansions not 
accurate enough and they are not generally used for the determination of a^(/i^). On the other 
hand this relation is useful for the study of the uncertainties which are introduced by the choice 
of particular scales; we will discuss this in the next section. 

Another alternative to the numerical evaluation of the coupling which is often used [19] is the 

„2 

expansion of the solutions in inverse powers of L = In so that the first term in Eqs. 1.4 to 1.6 
is the LO solution and higher terms are found recursively. Again, there is a certain freedom in the 
definition of A; up to NNLO we choose: 

2^_ 1 feihiL b\{\n^L-\nL-\) + b2 

where = ^ ■ As mentioned, although theoretically consistent, these approximate asymptotic 
solutions differ considerably from the exact ones at low scales (say, /^^ <2 GeV^). Since, as we 
will see, this region is relevant for the dynamical approach to parton distribution functions we 
will use the exact iterated solutions of Eqs. Eq. 1.4 and Eq. 1.6 all throughout this work. 



1.2 Factorization and the RG Evolution Equations 

The application of pQCD relies in what are known a.?, factorization theorems [6-10], which prove 
that factorization is possible for particular processes and observables, i.e., that the singularities 
which appear in the calculations either cancel or can be absorbed in the definitions of physical 
quantities. As we have seen, the ultraviolet divergences stemming from virtual radiative cor- 
rections are contained in the running coupling constant and are automatically included in the 
renormaUzed theory by changing the bare coupling for the miming coupling as{ix^). Depending 
on the gauge used in the calculations other ultraviolet divergences may emerge but they cancel. 

In general, they appear also infrared divergences with different (or mixed) origins: soft, related 
to the emission of gluons with vanishing energy — > Oj , or collinear/mass singularities, related 
to coUinear radiation off massless partons. If masses are introduced, the soft singularities remain 
while the coUinear ones are regularized given rise to a logarithmic dependence in the masses. 
After the combination of all (real and virtual) contributions, the soft divergences cancel while 
the remaining coUinear singularities are included in the definition of the parton distributions at 
a certain (arbitrary) finite mass scale known as the factorization scale jXp- As part of the proof 
of factorization theorems precise definitions of the parton distribution functions f (where /, 
stands for the distribution of the different kinds of partons, i.e., massless quarks and gluons) 
and coefficient functions (or partonic cross-sections) are formulated (for details see, for instance, 
[6-10]). 

A major point of factorization is that the infrared divergences appear at the partonic level in 
an universal way, i.e. with independence of the scattering process considered. Therefore the 
factorization scale dependence of the parton distributions is also of universal validity, it reflects 
the independence of physical quantities on fip and it is given by the evolution or renormalization 
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group equations (RGE): 

where the Bjorken-x is (related to) the longitudinal momentum fraction of the partons. The 

evolution kernels Pij or splitting functions stem from the colhnear divergences absorbed in the 
distributions and represent the (coUinear) resolution of a parton i in a parton j. They are calculated 
perturbatively as a series expansion in as{[Xp): 

oo 

Pij{§,i^j)^Pij{IMi^h)^Lcir\i^j)p^;\^-) (1.10) 

m=0 

At the moment they have been calculated up to 3 loops (m = 2), which is the necessary accuracy 
for a NNLO analysis; we summarize the expressions that we use in Appendix B. 

Since we work in fixed-order perturbation theory (m < 2), the last equality holds only up to 
higher order terms and the splitting functions (and therefore the PDFs) depend effectively on the 
renormalization scale as well. We can make explicit this dependence by expanding as{{Xp) in 
terms of as{ix^)', using Eq. 1.7 we get up to NNLO: 

=';f»>(»'l)+(Pi''-ft<'ln|)»?(Mi) 

+ ('',?'- (2ft4" + ft<')'"S + fi?''«''"^i)'''(''«)' 

which is useful to investigate the uncertainties due to variations on the renormalization and fac- 
torization scales separately. However, all available sets of parton distribution functions have been 
generated using iip = iir and we will also set them equal. 

The coefficient functions (also known as Wilson coefficients or simply as partonic cross- 
sections) are in general different for different processes but depend only on short-distance inter- 
actions and therefore are calculable in pQCD. They are expanded in analogy to Eqs. 1.10 and 1.11 
(with different starting powers depending on the process) and thus depend, besides on the phys- 
ical energy scale (/l^), on the factorization (/i^) and renormalization (/i|) scales as well. In an 
obviously symbolic notation, and writing only the dependence on the different scales, a general 
observable is calculated in pQCD as: 

C7(/i') = ^C,-(lng,a,(M|)) 0PDF,-(a,(M^)) 

In ^ , In g , fl,(M|)) ® PDF,- (in g , ^^(mI) ) 

where / denotes a particular contribution (subprocess) and PDF, represent the appropriate com- 
bination of PDFs to be combined with the partonic cross-section C,-. 

Beyond LO there is a certain freedom on where to include the finite contributions so that a 
factorization convention or factorization scheme is needed for a complete definition of the factors 
appearing in the factorization formulae. The most extended prescription is the so-called MS, 
where in addition to the divergent pieces only an ubiquitous ln4;r — Je is absorbed into the PDFs. 
Another commonly used scheme is the so-called DIS scheme, in which all the (light quark flavor) 
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finite contributions to the nucleon structure function F2 (Sec. 1.6) are absorbed into the parton 
distributions. We discuss the factorization scheme invariance and the transformations between 
these two schemes in Sec. 1.8. Global fits in both factorization schemes are presented in Sec. 2. 

By solving the RGE we determine the parton distributions at a certain scale from the distribu- 
tions at another scale and therefore we can make predictions for all 11^ once the x-dependence of 
the functions has been extracted from experiment at one scale known as the input scale. This is 
usually done by parametrizing the x-dependence of the parton distributions at the input scale and 
determining the parameters by means of global QCD fits, i.e., fits to significative data sets which 
(ideally) constrain all the PDFs for all x. 

An alternative to the direct numerical integration of the RGE Eq. 1 .9 is to work in the so- 
called Me//m n-space. The symbol (g) in Eqs. 1.9 and 1.12 denotes the Mellin convolution, which 
reduces to multipUcation for the Mellin («-)moments of the quantities involved, i. e., 

a"= C dx}^-^a{x)^{a®bf = aV. (1.13) 
io 

Hence, in the Mellin space, the evolution equations reduce to ordinary differential equations 
which are solved analytically; the x-space solution may also be obtained via a numerical Mellin 
inversion. On the other hand, this technique requires analytic continuations of non-trivial func- 
tions to complex values of n. We use the Mellin space to solve the RGE and for the calculation of 
nucleon structure functions. Other quantities, for which the moments of the coefficient functions 
are not known, are calculated in ordinary x-space. 



1.3 The Flavor Decomposition 



The calculations appearing in pQCD and in particular the RGE can be greatly simplified by 
expressing them in terms of particular combinations of parton distributions. For instance, the 
splitting functions are constrained by charge conjugation and flavor symmetry invariance, so that 
they can be written as: 



^ gm — 
P- — 



P 

^gq 
p 

^ig 
Inp 



' mi 



P_ _ = ^. pV I pS 



' mj 



= S P^- + P 



qq 



(1.14) 



where np is the number of quark flavors entering in the evolution (massless), which is not neces- 
sarily equal to the number considered in the miming coupling [20]. 

In view of the general structure of the splitting functions, it is possible to introduce some 
flavor combinations in order to decouple the RGE as far as possible. Note for instance^: 

dqr 



q =qi±qi^< 



dXnjJLp 

dqt 



pS np 

k=\ 



, where < 



np 



np 



pS 



MP^q- 



Pps=np{P^, + I^^) 



(1.15) 



^We work in the Mellin space although we do not explicitly write the n-dependence of the moments. 
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The minus combinations decouple from the gluon, whereas the combination maximally cou- 
pled to it: 



qs = q^ ^Y^qt =^ = Pqq^s+PqgS , where Pqq = +Pps 

k=i om^p 



(1.16) 



this combination together with the gluon are known as the (flavor) singlet parton distributions. 
Defining the vector q and a matrix of splitting functions P, they evolve as: 




I ^ - =Pq (1.17) 



Simple non-singlet combinations which evolve independently are: 

^NS,,7 = - ^ = Pm ^NS,o- ' with J, /: = 1 , . . .«f 

^NS = = L ^ = ^NS ^NS > where P^^^ = ^ns + ^NS 

k=i omfj,p 



(1.18) 



However, together with we will use for the evolution in the non-singlet sector even more 
convenient combinations given by [21]: 

k=l h^F ^j^g^ 

1 1 1 

or inversely: qf = ^ ,± - .v^ ^ + ^ L ^ J^^^t-, 

In the power expansion, starts at leading order while P^, P^^ and start at NLO, however 
I^q = P^q at this order, being different for the first time at NNLO. This allow for some simpli- 
fications at lower orders. As mentioned, at leading order Pqq^\ Fqq^^ and P^^^^ vanish and the 
structure of the splitting functions reduces to: 



pS(0) _ pS(0) ( \ p _pV _p± _ pV{0) ^ p(0) 



(1.20) 



At NLO P^J'^ =liq \ i-e-, Pqq=Pqq (hcncc Pns =^ ^NS = ^Ns) simplifications 
still hold. For instance, in view of Eq. 1.15 it is clear that each valence distribution (w^, d^) 
evolves independently and any initial symmetric sea (e.g., s~ = 0) is preserved by the evolution 
up to NLO. Beyond NLO we find the complete structure of the evolution equations where the q^ 
combinations evolve together in what is called the non-singlet sector and the q^ combinations 
evolve in conjunction with the singlet distributions in what is called the singlet sector. 
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1.4 Solutions of the RGE 

Since it is clear how to recover the explicit dependence on the renormalization and factorization 
scales separately, we will set here Hp = ;u| (henceforth simply ;U^) for simplicity in the discussion. 
With the definitions of last section all the evolution equations are of the same form, namely: 

^=^.(M^)^(M^) (1-21) 

where for the singlet distributions this equation is a matrix equation and one has to take care of 
the commutativity of the subsequent factors. 

Because all of the scale dependence of the splitting functions appears through as{{JL^), it is 
convenient to combine this equation with Eq. 1.2 and use as=as{n^) as independent variable. At 
leading order we can integrate immediately;formally: 

^ = -Ak^f^ q{as{n')) = {t)'''>o^L{^s,ao)qo, ^0,, ^ (1-22) 

where ao = as{lXg) and qo = q{jXo), being jx^ some reference scale at which the distributions are 
known. For the singlet distributions Ro is a matrix and the above solution is only symbolic. 
We can find an explicit solution by diagonalizing; the eigenvalues, projection matrices and some 
useful relations are: 

Eigenvalues: = ^ ^i?o,9? + Ro,gg ± \J {Ro,qq-Ro,ggy + 4Ro,qgRo,gq^ 

Projectors: = ^ => (R^f = R^, R++R- = l, R+R-=R-R+ = ^^'^^^ 

For an analytic function /: /(Ro) = /('"''") R"*" + /(''") R~ 

The last relation holds in particular for L{as,ao), which immediately gives an explicit solution for 
the singlet distributions at LO. 

To solve the equations at an arbitrary fixed order m [21, 22] we combine the evolution equa- 
tions at this order with Eq. 1.2 at the same order m: 

^ = -^fl<0^(«^)' Rm,q^'-^-tbiRm-i,q (1-24) 

where a.ga.iii bm = Next, we make a. series expansion around the LO solution as: 

oo 

q{as) = Ug{as)L{as,ao)U-\ao)qo, with Ug{a,) = a'^U^"'\ = 1 (1.25) 

m=0 

so that the problem is now reduced to the one of finding the operators By inserting this 

ansatz in the evolution equations we get: 

oo / oo \ / CO \ 

ufho,q+ I ar(f/f '^0,.-mC/f ') = /?o,,+ I a'^R^,, ^ d^ulT^ , (1.26) 

m=l \ m=l / \m=0 / 
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which order by order gives the commutation relations: 

m 

[t/f Uo,,] = Rm,q + m t/f ) , Rm,q = £ Rk,c, , (1-27) 

k=l 

which determine the Uq"^ operators. For non-singlet combinations the commutators vanish and 

(m) R 

the operators are directly Uq = — For the singlet combinations we can use the properties of 
the projection operators (second line in Eq. 1.23) to express them as: 

UW = -i(R+R„R+ + R R^R ) + ^^^"^ + ^"^"^^ (1 28) 

m r—r^—mr^ — r—m 

In addition to higher-order in the expansions, redundant terms (i^{a") with n>m) arise from 
the multiplications in the ansatz. One can remove these terms in advance by expanding in powers 
of Uo the factor U^^{ao) and carrying out the multiplications keeping only terms up to the desired 
order for consistency with the fixed-order perturbative approach. Up to NNLO we get: 

(1.29) 

where, again, for the singlet distributions one has to take care of the order of the different factors. 
This defines what are called the truncated solutions, which are the ones employed in our analyses 
since, as mentioned, they are more consistent with the fixed-order perturbative treatment. 

We have tested our singlet and non-singlet evolution codes using the PEGASUS program [23] 
for generating the truncated solutions together with the commonly used toy input of the Les 
Houches and HERA-LHC Workshops [24, 25]. For 10"^ <x<0.9 we achieved an agreement of 
up to four decimal places in most cases, which is similar to the required high-accuracy bench- 
marks advocated in [24, 25] . 

Note that the truncated solutions are not directly comparable with x-space evolutions or the 
so-called iterated n-space solutions, where all the terms are kept. Since the differences are 
of higher orders, whether or not to include which redundant terms is somehow a matter of taste. 
Notice, however, that these differences are not neghgible, specially at NNLO, where the precision 
aimed is of a few percent^. 



1.5 Isospin Violations in the Nucleon 

Electroweak radiative corrections may be included together with the QCD description of the 
structure of the nucleon [26-29] . In complete analogy with the emission of gluons, the emis- 
sion of photons by quarks introduce new mass singularities which are factorized in the parton 
content of the nucleon. The evolution equations of Eq. 1 .9 are consequently modified to include 
electromagnetic splitting terms; in general: 

^fl^ = I {Pij,QC^®fj + Pijm^®fj) (1.30) 

''For instance, using the toy input mentioned, we get for the NNLO (NLO) iterated/truncated ratios at ll^ = 10 GeV^ and 
;c=10-', 10"^ 0.1 respectively, the values 0.973(0.940), 0.992(0.975), 1.001 (1.002) for the sea quark {u + d) distributions and 
of 1.043 (1.049), 1.006 (1.016), 0.998 (0.996) for the gluon ones; the differences are smaller in the non-singlet sector. 
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where the sum over partons includes now additional distributions, e.g., the photon parton distri- 
bution function /(x,/^-^). 

The modified evolution equations contain additional mixing between distributions and flavor 
symmetry violations since different quark flavors have different charges and therefore evolve 
differently in what concerns the QED terms. In particular the standard isospin symmetry relations 
for the parton distributions of the nucleon, i.e., vanishing "majority" and "minority" asymmetries 
defined via 

dd{x,fx^) = dP{x,ix^) -u"{x,ix^), 

and similarly for the antiquarks 5q, are violated. Note that these asymmetries must observe the 
quark number and momentum conservation rules, e.g., separating, as usual, between valence and 
sea contributions, the first moments of the valence asymmetries must vanish. 

Although in principle one should solve the modified evolution equations, it is possible to get 
a first estimate, sufficient for most purposes, by working in LO QED; from Eq. 1 .30, the leading 
corrections are given by: 



d 



(1.32) 



7(0) 

where e/ is the electric charge and Pqq = similar evolution equations hold for the isospin 
asymmetries of sea quarks 5u and 5d. The evolution equation for the gluon distribution remains 
unaltered at this order. Notice that the addition [29, 30] of further terms proportional to ae^Pqy® y 
would actually amount to a subleading ff{a}) contribution since the photon distribution y{x^ii^) 
of the nucleon is of &{a) [31-34]. 

In principle one could expect the ^(a) corrections to be of the order of NNLO QCD; however 
the effects turn out [26, 28, 29] to be rather small (at most 1% in the relevant regions) and 
irrelevant for usual QCD applications. Furthermore they often mix with non-perturbative effects 
like differences in quark masses. Nevertheless they may be relevant for high precision experiment 
in which they appear somehow isolated. 

Since usually the data used in global QCD fits are rather insensitive to the isospin asym- 
metries, some plausible contour conditions for the above equations, i.e., some reasonable (non- 
perturbative) value for the asymmetries, must be assumed in order to integrate them. In this way it 
is possible to express the (LO QED) isospin asymmetries as functions of usual isospin symmetric 
parton distributions. We make such estimates in Sec. 2.1 1 within the framework of the dynamical 
model and in coimexion with the so-called "NuTeV anomaly", where they turn out to contribute 
significantly. 
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1.6 Unpolarized DIS Structure Functions 

Lepton-nucleon deep inelastic scattering (DIS) is the most powerful available tool for the de- 
termination of nucleon PDFs. DIS cross-sections are described in terms of standard kinematics 
[6-10]; we recall here only the definitions of the variables that we use (cf. Fig. 1.0): 

s={P + kf, W^={P + q)\ Q' = -q^, ^=S^q^ ^=9^' ^^'^^^ 

W is the invariant mass of the recoiling system, x is known as the Bjorken-x and y is often 
referred to as the inelasticity. For fixed center-of-mass energy (5), there are only two independent 
variables and the following relations hold: 

W^ = M^ + ^{\-x), Q^ = xys, (1.34) 

where M is the mass of the nucleon and the DIS regime implies ^ (deep) and » 
(inelastic). 

The differential cross-section for DIS can be written as products of leptoiuc and hadroiuc 
tensors. The leptoiuc tensors describe the lepton-gauge boson vertex and are well known [6-10]. 
The hadronic tensor for unpolarized DIS is usually parametrized as: 

W^y{P,q) = {-g^^ + ^)Fi{x,Q^) + ^F2{x,Q^) + ie^yap^F3{x,Q^) (1.35) 
q'^ F q ^ 2F q 

with Pn=P^ — -^q^. An alternative is to expand it using the polarization vectors of the (virtual) 
gauge boson; one gets then the so-called longitudinal structure Junction: 

Fl{x,Q') = {l+4-^x')F2{x,Q^)-2xF,{x,Q^)^F2{x,Q^)-2xFi{x,Q') (1.36) 

Q 

which is frequently used instead of T^i . A (sometimes) more convenient alternative definition of 
the structure functions is: 

^1,2,3 = (2Fi, -,^3) ^^L=- = ^2-^i (1.37) 

X X 
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(7+ + l^y" f )Fi - y^F[ ^ F-xFj) (1.38) 



The differential cross-section for DIS is given in terms of the structure functions as: 

dxdy xyQ^ 

where 7± = 1 ± (1 — y)^, / = NC,CC and — is taken for an incoming antiparticle (Z+ or v) and 
+ for an incoming particle (l~ or v). The factors stemming from the couplings and propagators 
of the gauge bosons have been included in the structure functions as ratios to those of the photon 
[35]. For charged currents: 

FL2M3=^^iiQ')FL2X3, with Kw {Q') = KwH , '^«'H = 2^ (1-39) 

where t] =2 for unpolarized charged lepton scattering and T] =4 for neutrinos. For NC: 
7NC _ z-r _„ ^^(n2\ rT^ 1 ^,2_l_^2^^2/n2^ ^--'^'^'^ - v--''^-^ — £ 



= -a, KziO") Ff + Iveae 4{Q^) 7f 



Kz{°°) 



GfmI 



-2 + 2sin^0iy) <2e = —\ 

(1.40) 

The Z° and contributions contain suppression factors which make the neutral-current cross- 
section being dominated by electromagnetic (7) exchange except at high Q^. Furthermore, the 
longitudinal structure function contributes sizably only for large inelasticity; a quantity which is 
often used is the "reduced" cross-section: 

a-. (?^y.)"'^=f,--gFrTfe.F3'«^ (1.41) 



dxdy 



which, as mentioned, approximately equals F^ for most of the data used in our analysis (Sec 2.2). 

In DIS there is a natural energy scale which by definition (deep) is much larger than the 
masses of the quarks flavors u, d and s (therefore called light), while it may be comparable to 
the masses of the c, b and t quarks (heavy). It is therefore convenient to consider separately the 
contributions to DIS structure functions stemming from light and heavy flavors. 

For the light quark contributions a massless description (in particular the massless evolution 
of Sec. 1.4) works fine and provides an appropriate natural choice for the factorization and 
renormalization scales; we give more details on their calculation in Sec. 1.7. The situation is not 
that clear for the heavy quark contributions since, besides Q^, the heavy quark masses appear as 
well as natural energy scales in the problem; we discuss these contributions in Sec. 1.9. 

Furthermore, since, as already mentioned, the DIS regime imphes (2^3>M^ and W^S>M^, the 
mass of the nucleon is often neglected. Its effects are only sizable at low Q^, say Q^< 100 GeV^ 
and medium to large x. In this region modifications to the usual (M=0) expressions or target mass 
corrections (TMC) must be taken in to account. The well-known expressions for the dominant 
"fight" F2^^ structure function are given in n-space by [36]: 

F2,TMc{n,Q)-^^[^) ^y(;^(„ + 2£)(n + 2£-l)+^lle^J ) ^^'^^^ 

where higher powers than (^)^ are negfigible for the relevant x < 0.8 region, as can straightfor- 
wardly be shown by comparing the above expression with the exact one in x-space [36]. 
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1.7 Light Quark Contributions to NC Structure Functions 

After factorization, the (leading-twist) light quark contributions to the neutral-current DIS struc- 
ture functions are given as^: 

where the sum runs over all the light quark and antiquark flavors, q'q wee the appropriate elec- 
troweak coupling factors and we have obviated the renormaUzation scale dependence. The coef- 
ficient functions are expanded in power series as usual; for j = 1 , 2, 3: 

oo oo 

C,-,..(a.) = E a^cj-) , Cfl = l; Cj,,,{a,) = aTc\rJ^, cfl = (1.44) 

m—O m=0 

where m coincides in this case with the number of loops entering in the calculations^. 

Combining the electroweak factors related with the exchanged boson and the quarks respec- 
tively, the charge factors appropriate for the complete NC structure functions are: 

4 = a'/'"- = el-2eqVeVqKz{Q^) + {vl + al)4{Q^), at = a+ 

=^1= -'^eqUq Kz{Q^) + AveaeVqUq K^{Q^) , = -a' , with: 

(1.45) 

= 5, V(7 = 2 — ^ sin^ 0(7 = 2 
ez) = — 5, = —j + J sin^ 9w, aD = 

where U(D) refers to u(d)-type quarks. Note that the couplings for quarks and antiquarks are 
either equal (1,2,L) or opposite (3), this allows for the rewriting of the structure functions as a 
sum over the quark flavors: 

np np 

!=1 !=1 

where np is the number of flavors entering in the evolution, i.e., light. To express them in terms 
of the flavors combinations introduced in Sec. 1 .3 we consider the following identity: 

np ± Up i np 

;=1 (=2 k=\ i=2 

Defining as well 0'*'= the appropriate combinations of couplings reproduce (upon nor- 

malization) those of the distributions in Eq. 1.19. 

It is important to note that the coefficient functions for quarks contain a pure singlet part 

stemming from interactions which cancel for the non-singlet combinations (Cj^^^=Cy ns for non- 
singlet), and a non-singlet part which contributes both for singlet and non-singlet combinations 
iCj^qq=Cj^s=Cj,NS+Cj^PS for singlet); redefining the gluon coefficient function as Cj^G^^npCj^qg 

^ We work in Mellin space although we do not explicitly write the «-dependence. Recall also that we have set /J^ = n| = 
*For the longitudinal structure function cj"^^ = c'^^^ - cf j^^ = 0; therefore Fl appears at the one-loop level, i.e., at LO is already 
of order Oj, which is in agreement with the naive parton model expectation F2=2xFi known as the Callan-Gross relation. 
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and using the above identity we get^: 

np 
i=2 

If 

= C;,NsK?5 + +a'^{Cj,Psqs + Cj,G8) (1.48) 

!=2 

np 

■^f^ = C3,NS {a'qls + L «i2-l^i2-l) 

The combinations stemming from one flavor only {qi) are given in this notation for the coefficient 
functions by: 

= 4 ifj^qt + 2C,,,,g) = a+ (C,,NS qj + ^-^qs + ^g) 
^3^0,?; ^ a-Cz,qqqi = a-.C3,NsqT 

At the moment all the necessary coefficient functions for a NNLO analysis are known: up to 
one-loop we use the well-known expressions, summarized for example in [37], and for the 2- 

(2) 

and 3-loop (for Fl) coefficient functions we use the parametrizations of [38, 39], except for C3 ,(,5 
for which we use the parametrizations given in [40]. 

1.8 The DIS Factorization Scheme 

The DIS factorization scheme is defined by the requirement that the light quark contributions to 
the function F2 retain their leading order form at all orders in perturbation theory. We can derive 
transformation relations between factorization schemes by comparing expressions for physical 
quantities, which must be, order by order, factorization scheme independent [41]. 

In order to manipulate simultaneously the contributions of quarks and gluons, it is convenient 
to regard the contributions to the structure function F2 stemming from a combination of quarks 
(q), and its related gluon contribution, as the quark sector of the product of a matrix of coefficient 
functions C2 times a general (quark and gluon sectors) distribution q; in analogy with the notation 
introduced for the singlet sector in Sec. 1.3. Since there is no corresponding gluon sector in the 
structure functions, only the upper part of the products is to be considered, i.e., beyond LO the 
gluon sector of the matrix of coefficients is undefined at this point. Keeping this in mind the 
mentioned contributions are given as: 

00 

^|°cC2q, C2=2K4'"\ Cf = l (1.50) 

m=0 

Combining this with Eqs. 1.2 and 1.9 and setting itt^=iU|= for simplicity, we can express its 
derivative as function of the structure functions itself: 

J^-C'S^c.-'+p+icPic-)^!, (1.51) 



'Actually the standard notation is somewhat more complicated, as for the splitting functions, there appear "+" and "-" combina- 
tions of coefficient functions. Since for NC only the "+" ("-") combinations are relevant for f2,L (^3) we omit these signs. 
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where P denotes the matrix of splitting functions appropriate for q. This makes clear that the 
quantity in brackets must be invariant under factorization convention transformations. Hence, if 
the coefficient functions are changed so that C2 +AC2 , the splitting functions must be 

conveniently changed as p('")^ AP^*"'. Expanding up to &{a^) and requiring invariance: 

ap(i) = m4'^-[ac(^\p('')] 

AP(2) = -2/3oC^^Uc5^' -^oAC^'^AC^'^ + 2Mcf +^1AC5^) -PWC^'^AC'^^ 

+AC(^Uc(^'p(0)-Ac(^)p(0)AC(^) + C^i'AC(^)p(0)-AC('^P(i) ^^'^^^ 
+p(i)A4^)-[ACf ,P(0)] 

On the other hand, the structure functions themselves must be (order by order) invariant un- 
der factorization convention transformations. Hence the parton distributions must change up to 
NNLO as q — > q + UgAq^^^ + a^Aq^'^); again expanding and requiring invariance: 

Aq(i) = -AC^'^q, Aq(2) = (C^'^AC^^^ +a41)aC^i) - ACf )q (1.53) 

As mentioned, the DIS factorization scheme is defined by C2= 1 at all orders, hence the 
transformation of the F2 coefficient functions from any other scheme (in practice from MS) is 
given by AC2^^ = c'^^^^^^ _ c^™^ = — Cj'"'' , Vm > 1 . The general results for the splitting functions 
and the parton distributions reduces then to: 

AP(1) = -A)C?) + [C?\P(0)] 

AP(2) = ft4')cW - 2A)C^'' - ^C^'^ - [4'\P(0)]C^'' + [C^'',P(i)] + [Cf ,P(0)] (1.54) 
Aq(i) = c(i'q, Aq(2) = c(')q 

These relations refer to a general distribution q, it is convenient to make explicit the trans- 
formations for the flavor combinations introduced in Sec. 1.3. For non-singlet combinations the 
gluoiuc contributions cancel and therefore the above relations may be applied directly to get the 
respective non-singlet splitting functions and distributions in the DIS factorization scheme; note 

that the commutators vanish. In the singlet sector, however, these relations do not completely fix 
the transformations because the gluon sector of the matrix of coefficient functions is still undeter- 
mined. The conventional choice is to extend the restriction on the coefficient functions stemming 
from momentum conservation to all moments; this results in: 

, , / dl \ / d'J d"i \ 

\ '-im ^"-P^iml \ 2,5 ^ifil 

(1.55) 

Since in order to invert the transformations a change ACj = is needed in the coefficient 

functions, the same expressions with the appropriate sign changes give the inverse transformation. 

1.9 Heavy Quark Contributions and their Resummation 

The contributions to hard scattering processes involving heavy quarks h = c,b,t become increas- 
ingly important above the threshold for the production of a quark-antiquark pair of heavy flavor 
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mnnnf 




Figure 1.2: ff{as) (leftmost) and some illustrative ff{(xj) diagrams contributing to the virtual photon-gluon 
fusion process y*g hh. Note that beyond LO other production mechanisms, e.g. Y*q hhq, etc. appear The 
figure has been taken from [79]. 



(e.g. > 4m^ in DIS, where m stands for the mass of the heavy quark). Although in principle 
there could exist an intrinsic (initial-state) heavy quark content of the nucleon, all data are well 
described through extrinsic heavy quark production only (Sec. 2), i.e., being all the heavy quarks 
generated from the (initial-state) hght quarks {u, d, s) and gluons, so that it appears that the in- 
trinsic heavy quark content of the nucleon, if any, is marginal. The extrinsic production of heavy 
quarks is in principle calculable in fixed-order perturbation theory and is known as the fixed fiavor 
number scheme (FFNS) since, in contrast to variable fiavor number schemes (VFNS) to be dis- 
cussed below, the number of light quark flavors (partons) does not change with the factorization 
scale. 

In the case of the neutral-current (we drop here the superscript NC) DIS structure functions 
of Sec. 1.6, the (subleading) heavy quark contributions enter as: 



plight I pheavy 



2.L 



riheavy 
2,L 



b 

2,L 



(1.56) 



where "light" refers to the contributions of Sec. 1.7 and top contributions are negligible at present 
energies. Their contributions to F3 vanish at LO [42] and are negligibly small at higher orders^. 

The ^(ots) FFNS contributions due to the virtual photon-gluon fusion subprocess y*g — > hh 
are well known [42-45] and their QCD corrections (cf. Fig. 1.1) have been calculated so far up to 
ff{ccj) [46, 47]. Although complete analytic expressions are not available for all the coefficient 
functions, the results are contained in computer programs. In particular we use the code offered 
in [47], which combines some known analytic expressions together with grids for the more com- 
plicated coefficient functions; they are represented in the MS factorization scheme. Up to NLO 
the relevant heavy quark contributions to DIS structure functions are given by^ : 



F,l2,Li^,Q^m^j-- 



An^m^ 



z 



{atcfl[^,^)g{l,n') 



.(1) 
-Kg 



(1.57) 



+a+(4;j(r7,^)+4'^)(r7,^)lng)^(f,Ai^))]}, 



^This can be guessed by observing that at the relevant large values « /; — h, where the meaning of the ejfective heavy quark 
distributions h will be clarified below. 

'We use the expressions in [47], which were slightly modified respect to those in [46] due to an additional mass divergence which 
appears when the virtual photon coupled to the light quark goes on-mass shell, leading to the appearance of the function 4',' (4',' =0)- 
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where the sum runs over all the light quark and antiquark flavors and the charge factors 
appropriate for the complete NC structure functions are given in Eq. 1 .45. Note that the coefficient 
functions indicated by c's (d's) originate from partoiuc subprocesses where the virtual photon 
is coupled to the heavy (hght) quark and that those indicated by a bar appear through mass 
factorization [46, 47]. The scaling variables (t],^) are defined as: 




(1.58) 



being s= ^(1 — z) the square of the cm. energy of the virtual photon-parton subprocess (i.e., 
the partonic version of VK^, cf. Eq. 1.34). With these definitions the threshold region (s ~ 4m^) 
is characterized by t] ^ 0, while the asymptotic regime (Q^ S> m^) is given as ^ ^ oo. Notice 
that these are two independent limits. In principle one can consider the asymptotic structure 
function defined by F^'''^^™^(x, 2^,m^) = limQ2^^2F^{x,Q^,m^), however, in contrast to what 
is sometimes believed, even in this limit F^^ receives contributions from the threshold region as is 
clear from Eq. 1.57, where the collinear momentum fraction (z) is integrated out, i.e., s takes all 
its possible values 4m^ <^<x(^^'^) — particular those close to the threshold. 

Further, /i in Eq. 1 .57 stands for the factorization scale, which has been set equal to the renor- 
malization scale for simpUcity (recall Sec. 1.2). A dedicated study [48] of the perturbative sta- 
bility of heavy quark production at high energy colliders concluded that the factorization scale 
in Eq. 1.57 should preferably be chosen to be jX^ = Am^ at LO, while the NLO results are rather 
insensitive to this choice (usually either /i^ = 4m^ or /i^ = 2^ + 4m^); even choosing a very 
large scale like /^^ = 4(2^ + 4m^) leaves the NLO results essentially unchanged, in particular at 
small-;c. Furthermore, they concluded that the fixed-order calculation for heavy quark produc- 
tion is entirely reUable and perturbatively stable, provided that one employ consistently parton 
distributions and strong coupling constant values of the appropriate order [48]. As has already 
been mentioned, the extrinsic heavy quark production is, in addition, experimentally required, in 
particular near the threshold region. 

The ff{oLf) 3-loop corrections to F/' and first rudimentary contributions to F2 have been 
calculated recently [49-51] in the limit Q^^np-. However these asymptotic results are neither 
applicable for our present investigations nor relevant for the majority of presently available data 
at lower values of Q^. There exist also soft gluon resummations [52, 53] which include the 
logarithmically enhanced terms near threshold up to NLL, which improve the convergence of 
the perturbative series in the very small-x region and gives a first approximation (NLO + NLL) 
towards the NNLO results. Our ignorance^^ of the full fixed-order corrections to 

constitute the major drawback for any NNLO analysis of DIS data^ ^ . 

The coefficient functions appearing in Eq. 1.57 contain terms with factors of the form In^ 
which in principle could vitiate the stability of the calculation and might suggest that these con- 

'"it should be mentioned that we have attempted to mimic the NNLO contributions in the relevant kinematical regions by naively 
assuming them to be down by one power of respect to the NLO ones and being given by a constant S'-factor times different 
combinations of the lower order expressions evaluated using parton distributions of different orders. The results (Sec. 2) were however 
insensitive to such ad hoc corrections and, furthermore, this approach (guess) appears to be not appropriate since playing the same 
game at NLO can not reproduce the correct results in the relevant kinematical regions. 

''Note however that for totally inclusive (light + heavy) structure functions the heavy quark contributions, although important, 
enter at subleading levels. 
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tributions should be resummed [54, 55]. In its simplest version, what is known as the zero mass 
(ZM) variable flavor number scheme (VFNS), the asymptotic FFNS heavy quark coefficient func- 
tions are separated into the massless Wilson coefficients of the light partons and massive oper- 
ator matrix elements which are used as transition functions to define new parton distributions 
for "heavy" quarks, which are subsequently treated as massless partons within the nucleon; i.e. 
whose contributions to DIS structure fimctions follows Sec. 1.7 and whose distributions h{x^ii^) 
are treated as those of the light quarks in Sec. 1 .4, being usually (up to NLO and choosing m = n) 
generated from the boundary conditions h{x,m^) = h{x,m^) =0. 

In general the matching conditions are fixed by continuity relations [54, 55] at the unphysical 
threshold }Xp = m^, where the "heavy" >3 quark distributions are generated from the np — I 
ones via the massless renormalization group equations of Sees. 1.2 to 1.4. Hence, this factoriza- 
tion scheme is characterized by increasing the number of massless partons np by one unit at the 
unphysical "thresholds" ii^ = iip = m^, in a similar way as is done in the MS renormalization 
scheme for the strong miming coupUng (Sec. 1.1); note however that, as already mentioned, in 
general np^fiR [20]. 

Despite its theoretical problems (e.g. the transmutation of a final-state quark into an initial- 
state one), a major general advantage of the ZM-VFNS is that it simplifies considerably the 
calculations, which in many situations become unduly complicated in the FFNS. For example, the 
single top production process at hadron colliders via W'^-gluon fusion requires in the FFNS the 
calculation of the subprocess ug ^ dtb at LO and of ug^ dtbg, etc. at NLO while in the ZM- 
VFNS one needs merely ub ^ dt at LO and ub ^ dtg, etc. at the NLO of perturbation theory 
[56]. Obviously this approach is convenient because of its simplicity, as is clear in particular 
for processes for which the FFNS coefficient functions containing explicitly the complete mass 
dependence are not known. It should nevertheless be regarded as an effective treatment of heavy 
quarks, keeping always in mind that the existence of an intrinsic heavy quark content of the 
nucleon is experimentally disfavored and that, being generated from the asymptotic heavy quark 
contributions, it cannot reproduce the fully massive result (FFNS), not even at very high [48] . 

There exist also more involved schemes with a variable number of active flavors (for a recent 
review see [57]), sometimes referred as general mass (GM) VFNS, where mass-dependent cor- 
rections are maintained in the hard cross-sections and a model-dependent interpolation between 
the ZM-VFNS (for the asymptotic regime) and the FFNS (near threshold as required experimen- 
tally) is achieved. Since the threshold behaviour of the FFNS is included in these models, they 
are generally able to reach a good agreement with data; in fact there are no observable signatures 
which allow to uniquely distinguish between these GM-VFNS schemes and the FFNS. On the 
other hand, the universality of the distributions is somehow put in danger due to the inclusion 
of non-coUinear terms in the partonic picture. Furthermore, the partonic cross-sections for most 
processes are either calculated in the simplest ZM-VFNS or already known for the FFNS, so that 
ultimately the user must choose one of these two schemes. Needless to say that the combination 
of GM-VFNS distributions with these expressions, despite being frequent, is inconsistent and 
should be avoided. 

In our opinion the stability of the FFNS renders attempts to resum supposedly "large loga- 
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rithms" in heavy quark production cross-sections superfluous. Since besides no intrinsic heavy 
quark content of the nucleon is needed experimentally, only the light u,d,s quark flavors and 
gluons constitute the "intrinsic" genuine partons of the nucleon and the heavy c, b,t quark flavors 
should not be included in its parton structure. For applications for which the FFNS expressions 
are not known, it is possible to generate ZM-VFNS distributions (based on FFNS fits) to be 
used consistently in these calculations at large enough scales. We compare FFNS and ZM-VFNS 
predictions for several illustrative processes in Sees. 2.9 and 2.10 (cf. [4]), and show that the 
ZM-VFNS can be employed for calculating processes where the invariant mass of the produced 
system is sizeably larger than the mass of the participating heavy quark flavor, for instance for 
LHC phenomenology. 

1.10 Hadroproduction of Vector Bosons and Jets 

In spite of its title, the aim of this section is not to provide a detailed discussions of hadroiuc 
cross-sections in QCD, moreover since the major experimental input for the determination of 
our dynamical distributions in Sec. 2 comes from deep inelastic scattering experiments. For 
completeness and reference we include here a brief overview of the hadron-hadron scattering 
processes that we use in addition to the DIS structure functions, namely Drell-Yan dimuon pro- 
duction as well as and inclusive jet production. 

As discussed in Sec. 1.2, a fundamental property of the QCD partonic description of the 
hadroiuc structure is that the colHnear divergences appearing in the calculation of any hard- 
scattering process (which originate the evolution equations) are universal, i.e., the same ones 
in all processes. Therefore the parton distributions appearing, for example, in the expressions 
for DIS structure functions are the same that describe the structure of the incoming hadrons in 
hadroiuc production. In the latter case the general factorization formula of Eq. 1.12 reduces to: 

o{Pi,P2)=Y^ / dxidx2dij{xi,X2M^,as{lx^) \ fi{xu^'^) fj{x2,^^) (1.59) 
hi 

where M is the relevant energy scale (e.g. the invariant mass of the muon pair for Drell-Yan, 
the mass, etc.) and ii is the factorization and renormalization scale, which have been set 
equal for simpUcity (recall Eq. 1.7). Again Cij is calculable in fixed-order perturbation theory as 
a series expansion in as{jJ?') which starts which different powers depending on the process. 

As in the case of DIS, hadron-hadron collisions are described in terms of standard kinematics 
[6-8]; we recall here some definitions (cf. Fig. 1.2): 

S={Pl+P2f, M^=S={xiPi+X2P2r ^XlX2S, X = — ^ X1X2 (1.60) 

the four-momentum of the produced system is conveniently given as: 

{E,px,Py,Pz) = (niT coshy, pt sin(j) , pt COS (p ,mT sinhy) (1-61) 
where pr is the transverse momentum, m\=M^ + pj and <j) is the azimuthal angle. The rapidity 
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Figure 1.3: Parton model picture of hadron-hadron scattering. The filled circles represent the structure of the 
proton and the rectangle the hard partonic interaction. 



y and, equivalently, the Feynman-x are defined by: 

which lead to the relations: 

;CF = 2(^T+4)'sinhy^y = arcsinh^i;cF(^T+4) (1-63) 

Depending on the detected final state it is possible to define different cross-sections for 
hadron-hadron collisions. In the Drell-Yan mechanism [58] a quark from one hadron and an anti- 
quark from the other one annihilate into an intermediate vector boson (7, or W^) which subse- 
quently decays into a lepton pair. Considering lepton pairs of invariant mass M<^Mz the process 
is dominated by virtual photon exchange. The experiments usually consist on proton (beam)- 
nucleon (proton or some other nuclei) collisions and the cross-sections are extracted from the 
detection of muon pairs (dimuon production) produced in the decay of the virtual photons^'^. 

The relevant LO/NLO double-differential distributions -j^i;. for the Drell-Yan pro- 

cess have been summarized in the Appendix of [59]^^. More recently the NNLO corrections to 
the rapidity distribution have been calculated as well [62, 63]. Combinations of parton distribu- 
tions like^"^ q\q2 + qiqi, etc. appear in these expressions and therefore the Drell-Yan process is 
sensitive separately to sea and valence distributions, which is in contrast to the NC DIS cross- 
sections, where they usually enter as ^ -I- ^ (cf. Sec. 1 .7). This is the main reason for the inclusion 
of Drell-Yan data in global QCD analyses of PDFs; in particular they are instrumental in fixing 
d/ u (or d — u), for example ~ ^ (l + ^) for xi S> X2 at LO. Another alternative is the use of 
CC (neutrino) deep inelastic scattering structure functions. 

A process closely related to Drell-Yan dilepton production is the hadronic production of elec- 
troweak bosons [7-9]. Here we briefly outline the calculation of total rates for production, 
in particular their dependence on the factorization scale, because we use them in Sec. 2.10 in 
connexion with the perturbative stability of different treatments of heavy quark masses at high- 
energy colliders. Since the decay widths of the electroweak bosons are small compared to their 

'^There are also prompt/direct photon experiments which detect a real photon with high pj; we do not consider them any further. 
^'There is an error in Eq. (A.8) of [59], which has to be modified [60, 61] in order to conform with the usual MS convention for 



the number of gluon polarization states, 2(1 — e) in 4 — 2e dimensions. 

'■ dMdxF 



'''To keep the notation compact we use here qi=fg{xi,n^) with ( = 1,2 and below in this paragraph 0?"= ^^J^ . 



25 



1. Elements of QCD 



masses, the narrow width approximation may be used to integrate the differential mass distribu- 
tion ^ [64]. Using an arbitrary factorization scale jX^ = 0'{M^) we get at LO^^: 

a^\s) = ^V2G^Tw Vl,^ C ^-^q,ixul^')qji^,ll') (1.64) 

where Tw = Vg.qj are the CKM matrix elements and the sum runs over all light quark and 
antiquark flavors in both hadrons, i.e. the contributions considered originate from partonic sub- 
process involving light flavors (ud^ W^, us ^ W^, etc.). 

The heavy quark flavor contributions to the total production rate in the ZM-VFNS are 
calculated in analogy to the light contributions (cs{d)^W^, bu^W^, etc.) while in the FFNS 
they are, like in DIS, gluon induced (i.e. proceed via gs{d)^cW'^, gu^bW^, etc.), for exam- 
ple, for the subprocess gs—>-cW^ they are given by: 

(7^^^"^(^) = f' dxi f' dx2 [g{xul^Mx2,l^^) + (1 ^ 2)] a^-'-^'^JxiX2s) (1.65) 

where now T= ("'c+^w) ^gs-^w found in [4, 65]. Unfortunately, the NLO corrections 

to this (massive) FFNS cross-section are not available in the Uterature. Only quantitative LO and 
NLO results for the analogous process gb^tW~ have been presented in [66], but questioned in 
[67]. The NLO/LO K- factor is expected [67] to be in the range of 1.2 - 1.3 . 

Another important process in hadron-hadron collisions is the production of high-/??- jets re- 
sulting from the hadronization of partons produced in hard scattering processes. Although non 
straightforwardly, it is possible to define infrared-safe cross-sections and to compare data and 
theoretical predictions [6-8]. Of particular interest for us is the pr dependence of these cross- 
sections, which fall strongly as pr increases. These distributions help to constrain the strong 
coupling and are specially important for the determination of the gluon distribution at high x, 
where no other experimental information directly sensitive to the gluon distribution is available. 

Since the detected objects are not directly the (colored) partons but rather the result of their 
fragmentation into bunches of color neutral particles (jets), any description of jet cross-sections 
is strongly dependent on the (Monte Carlo) algorithm used to define the jets; needless to say that 
it is crucial for an appropriate comparison that both data and theory use the same jet definitions. 
We use the so-called cone algorithm as implemented up to NLO in the f astNLO package [68], 
in which the relevant LO/NLO quantities needed for the calculation of jet cross-sections are 
tabulated for each particular data set in a way that permits the evaluation of the cross-sections 
using arbitrary parton distributions. Jet cross-sections beyond NLO have not yet been calculated. 



^In Sec. 2.10 we evaluate these rates up to NLO, however, since the expressions are somehow lengthy and cumbersome and, 
moreover, the corresponding IJpNS heavy quark contributions are only known at LO, we limit the discussion here to LO; which is 
sufficient for the exhibition of the main features. The reader interested in the NLO expressions is referred to [64]. 
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2.1 The Dynamical Parton Model 

The parton distributions of the nucleon are extracted from experimental data by two essentially 
different approaches^ which differ in their choice of the parametrizations of the input distributions 
at some low input scale /Xq. In the common approach, hereafter referred to as ''standard", e.g. 
[71-109], the input scale is fixed at some arbitrarily chosen I^-q > I GeV^ and the corresponding 
input distributions are unrestricted, allowing even for negative gluon distributions in the small 
Bjorken-x region [84, 85, 107, 108, 1 10], i.e., negative cross-sections like Fl(x, Q^). 

Alternatively [3-5, 111-119] the parton distributions at /i^>l GeV^ are QCD radiatively 
generated from valence-Mke, i.e. positive definite (xf >0 with xf ^ for x 0), input distri- 
butions for all partons at an optimally determined input scale Mo < ^ GeV^. This more restrictive 
ansatz imphes, as we will see, more predictive power and less uncertainties concerning the be- 
havior of the parton distributions in the small-x region at li^ > Hq, which is to a large extent due 
to QCD dynamics. 

The dynamical description of the nucleon connects the naive quark model with the parton de- 
scription of deep inelastic phenomena, i.e., the idea that baryons are bound states of three quarks 
with scaling violations and the fact that about half of the momentum of the nucleon is carried by 
the gluons. The original dynamical assumption [111-114] was that at a certain sufficiently-low 
resolution scale the nucleon consist of its constituent/valence quarks only, appearing the gluon 
and sea quarks as a result of bremmstrahlung processes, therefore radiatively or dynamically 
generated. 

In order to improve the agreement with constraints imposed by data on direct-photon produc- 
tion and some deep inelastic data at medium x (see [116, 1 17] for more details) the model was 
subsequently extended to include gluon [116] and sea {U = d) [117] valence-Uke input distribu- 
tions, amounting the underlying physical picture to "constituent" gluons and sea quarks comoving 
with the valence quarks at the low input scale Hq. The outcome of these modifications are steep 
gluon and sea distributions at small x which, as mentioned, due to the valence-like structure of 
the input distributions and the low value of the input scale (longer evolution distance), are mainly 
attributable to QCD dynamics. This is in contrast to the fine tuning (or even ad hoc extrapolations 
into unmeasured regions) required in standard fits, which depend strongly on the assumed ansatz. 

^ There are also early attempts [69, 70] to extract them using neural networks, which we do not consider any further. 
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These rather unique predictions were first confirmed at HERA [120, 121] (see also [122]), 
which appear to support the dynamical approach to the determination of parton distribution func- 
tions, i.e., that the low-x structure of parton distributions can be generated dynamically, almost 
in a parameter-free way, starting from valence-like distributions at a low input scale. In any 
case, it does not have any disadvantage^ respect to "standard" approaches, while its predictive 
power is an important and desirable feature. We will analyze quantitatively the implications of 
the dynamical approach to parton distribution functions in the next sections. 

Besides other criticisms [118], whether or not the standard evolution equations are applicable 
down to such small scales (;Uq ~ 0.5 GeV^) may be a matter of controversy. It should be noted 
that the predictions of the model are only intended to be compared with physical observables at 
higher scales, say /i^ > 2 - 3 GeV^; below that, higher-twist effects (although decoupled from 
the evolution of covariant leading-twist distributions themselves) may become relevant and even 
dominant. Nevertheless the dynamical model provides also a natural link coimecting nonpertur- 
bative models valid at < I GeV^ with the measured distributions at > 2 - 3 GeV^. 

2.2 Experimental Input and Analysis Formalism 

Following the radiative approach we have updated the previous LO/NLO GRV98 dynamical par- 
ton distribution functions of [1 19] and extended these analyses to the NNLO of perturbative QCD. 
In addition we have made a series of "standard" fits in (for the rest) exactly the same conditions 
as their dynamical counterparts. This allows us to compare the features of both approaches and 
to test the the dependence in model assumptions. Beyond LO we use the modified minimal sub- 
straction (MS) factorization and renormalization schemes. In order to test the dependence of our 
results on the specific choice of factorization scheme we have repeated our NLO dynamical anal- 
ysis in the DIS factorization scheme of Sec. 1.8. The quantitative difference between the MS and 
DIS results turn out to be rather small. 

The statistically most significant data that we use are the HERA (HI and ZEUS) measure- 
ments [106, 123-126] of the DIS "reduced" cross-section of Eq. 1.41^ for Q^>2 GeV^. Since 
the experimental extraction of the usual (one-photon exchange) fJ from is (parton) model 
dependent, we have chosen to work with the full NC framework in order to avoid any further 
dependence on model assumption. However, it turned out that fitting just to fJ gives very similar 
results. In addition, we have used the fixed target F/ data of SLAC [127] , BCDMS [128] , E665 
[129] and NMC [130], and the structure function ratios F^/F{ of BCDMS [131], E665 [132] and 
NMC [133]; both subject to the standard cuts Q^>4 GeV^ and > 10 GeV^. The well-known 
standard target mass corrections (Eq. 1 .42) to F2 have been taken into account in the medium to 
large ;ic-region for < 100 GeV^. 

^Except maybe a marginally larger for the statistically most significant data sets included in the fits due to the more restrictive 
ansatz; for the same reason the model is expected to fit comparably (or better) data with lower significance. 
'Note that the data we use are the radiatively corrected ones as presented by the experimentalists. 
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As indicated in Sec. 1.9, we work in the (three flavors) FFNS, i.e., heavy quarks (c.,b,t) are 
not considered as partons and the number of active flavors appearing in the splitting functions 
and the coefficient functions is fixed to nf=3. An extension to the usual (ZM)VFNS will be 
presented in Sec. 2.9. Note that in all our analyses the strong running coupling is governed by a 
variable number of flavors jir, as is common in the MS renormalization scheme. This is perfectly 
compatible with the chosen value for nf as detailed in [20]. 

The heavy quark (c, b; top quark contributions are negligible) contributions to f^^^ are theoret- 
ically described in the FFNS by the fully predictive fixed-order perturbation theory (cf. Sec. 1.9). 

yZ Z 

These contributions are quantitatively negligible for FjY' in Eq. 1.40, but they have been in- 

yZ Z 

eluded nevertheless; for ' they are negligibly small (cf. Sec. 1.9) and we have neglected 
them. The HERA measurements of heavy quark production [F2^) of [134-137] have also been 
included in the LO and NLO fits. In the NNLO fits these data have not been included, and the 
unknown third order coefficient functions in the heavy quark contributions of Eq. 1.56 have been 
neglected. For the quark masses we have chosen: 

Mc = 1.3GeV, m;, = 4.2GeV, m, = 175GeV, (2.1) 

which turn out to be the optimal choices for all our subsequent analyses, in particular for heavy 
quark production. 

Furthermore, the Drell-Yan dimuon pair production data of the E866/NuSea (fixed target) 
experiment for with N= p, d of [138] as well as their asymmetry measurements [139] 
have been used. The description of Drell-Yan cross-sections at LO of QCD is well known [6-8] 
to be systematically out of the ballpark by about 30 to 50%. This is taken into account in our LO 
fits by including a phenomenological (constant) K-factor (i.e., c^^^Ko^^) which is allowed to 
float, in both cases (dynamical and "standard") we obtain K=IA6 . A further complication in the 
inclusion of these data is that they are given in terms of xp, whereas the NNLO expressions have 
been given only in terms of the dilepton rapidity y [62, 63]. Since experimentally the dilepton 
Pt is small (below about 1 .5 GeV) as compared to the dilepton invariant mass M >5 GeV, we 
have checked that it can be safely neglected and the two distributions can be related using leading 
order kinematics^, as was done in [105]. For our NNLO analysis we used the routine developed 
in [63] and improved in [105]. 

Finally, the Tevatron high-/?! inclusive jet data of DO [140] and CDF [141] have been used 
together with the f astNLO package of [68] for calculating the relevant cross-sections up to NLO. 
In the NNLO fits the jet data have been omitted due to the lack of its theoretical description at 
this order. 

''That is, considering pj- = so that 3^ = ^^mTy (cf- Sec 1.10). 
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Figure 2.1: Kinematic x/Q^ range of fixed target and collider experiments; also shown are the major constrains 
that they make on the PDFs. The figure has been taken from [19]. 



The kinematical range of fixed target and collider experiments are complementary as can be 
seen in Fig. 2.0. Note however that in the theoretical description of the data usually convolutions 
of the parton distributions are involved, i.e., this range is only orientative of the values at which 
the parton distributions need to be evaluated. 

It is worth to mention that in the (NLO) DIS analysis the description of heavy quark contri- 
butions, Drell-Yan data and jet data is achieved, using their known theoretical MS expression, 
through the factorization scheme transformations of Sec. 1.8. These transformations allow also 
for a consistent comparison of our DIS results with the ones obtained in the MS scheme. 

We include as{M^) as a free parameter in our fits, and determine its value together with the 
parton distributions. It should be noted that there is a certain correlation between the value of 
the input scale jx^ and the resulting values for as{M^), which tends to increase with jx^. In the 
dynamical approach we performed fits for various values of the input scale keeping as{M^) as 
a free parameter while requiring a valence-like structure {a/ > 0) for all the input distributions. 
Then we fixed the best choice for ^u^ = 0.30, 0.50, 0.55 GeV^ for the LO, NLO and NNLO fits 
respectively and performed the final precision fits and error analyses. In the "standard" approach 
the input scale was fixed to ^Uq =2 GeV^. 

On the line of [118, 119] we use the following parametrizations for our [3] input distributions: 

XM„ {x,l^l^)=Nu x"" ( 1 - x) ^" ( 1 + A„ + x) 
xd,{x, IX^) =Ndx"'>{\- x)''^' (1 + + Bdx) 

x^{x,^l)=N^x"'^{\-x)^^{\+A^^ + Bl^x) (2.2) 
xI(x,Aio^) = A^zx°^(l -x)*5;(i +Aev^ + Bzx) 
xg{x,^l)=Ngx''>i{\-xfg 
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where = u — u, = d — d,l, = u + d and A = d — u. The distributions are further constrained 
by quark number and momentum conservation sum rules: 



which, as usual, we use to determine A^«,d,g, which therefore are not free parameters in our fits. 
We have tried different forms for these parametrizations, including (l+cx^) and (l+Cy/x+dx + 
ex^'^) for the polynomials, without finding any improvement or substantial change relative to the 
classical form Eq. 2.2. In particular all of our fits did not require the polynomial for the gluon 
distribution. The values obtained for the parameters of the input distributions in our different fits 
are given in Tab. A. 1 . 

As in [118, 119], the "light" sea is no longer symmetric {u^d) as required by the Drell- 
Yan data, which are instrumental in fixing A (cf. Sec. 1.10). Since the data sets we are using 
are insensitive to the specific choice of the strange quark distributions, we consider a symmet- 
ric strange sea. In the dynamical approach the strange densities are entirely radiatively gener- 
ated starting from 5(x,/io) = ■^(•^^Mo) ~ ^' while in the "standard" approach we choose as usual 
s{x, III) = s[x, Aio ) = ? ("(^, j"o ) + d{x, III)). 

The minimization procedure follows the usual chi-squared method with defined as: 



where p denotes the set of 21 independent parameters in the fit, including aj(M|), and is 
the number of data points included; AT = 1739 (1567) for the LO/NLO (NNLO) fits. The errors 
include systematic and statistical uncertainties, being the total experimental error evaluated in 
quadrature. The (fully correlated) normaUzation error of each data set is treated separately by 
allowing a common normaUzation factor (A^set) to float within the experimental error^. The values 
obtained for these normalization factors in our different fits are given in Tab. A.l. 

The minimum obtained in our fits as well as the contributions stemming only from subsets 
of data are given in Table A.l. As expected is slightly smaller for the "standard" fits, in 
particular for the DIS data which are the statistically most significant sets. In general, the LO fits 
do not give an appropriate description of the data while at NLO it is already satisfactory {x^ — 1) 
and the NNLO fits result in an even better (smaller) x^, typically xlmLo — 0-9Xnlo- compare 
extensively our results with the data used in Sec. 2.6. 

In view of this improvement in x^^ one could in principle reduce the number of parameters 
at NNLO and still obtain a good fit to the data {x^ — 1), which would result in a reduction of 
the NNLO uncertainties. As explained in the next section, we have chosen for the uncertainty 
estimations the (relatively simple) hessian method, keeping the same ansatz in all of our fits 

'These normalization factors multiply the data in our current analyses, i.e. data(!) — > iVset(i) data(/) in Eq. 2.4. In [3] a slightiy 
different definition of in which these normalization factors entered dividing the theoretical predictions, was used and thus the 
values quoted for in [3] are slightiy different from those in Tab. A.l. These difference are in any case negligible and we have 
explicitiy checked that this does not affect in any case the outcome of the fits. 




(2.3) 




(2.4) 
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in order to make these estimations somehow comparable. Note that the error calculations for 
parton distributions are rather arbitrary and must be interpreted with care and understood always 
as a rough estimation, useful only to compare analyses with similar treatments and in particular 
to disentangle in which regions the distributions are more/less constrained. This would be the 
case even with a more rigorous error analysis (taking into account correlations, etc.) since, as 
mentioned, the results obtained depend strongly on the framework used. 

Furthermore, these error estimations constitute the propagation of the experimental errors into 
the input parameters of the distributions, other sources of errors like, for example, the choices of 
renormaUzation and factorization scales, which particular solutions of the RGB are employed 
(cf. Sec. 1.4) or even which data sets are used, should be treated separately and, in any case, 
keep always in mind that they are not included in what is usually meant by "error" of the parton 
distributions. 

2.3 Estimation of Uncertainties 

The experimental errors of the data induce an uncertainty in the determination of parton distribu- 
tion functions (and therefore to any quantity calculated using them) which can be estimated now 
by linear propagation and other methods [142-146]. 

Our estimates are based in the so-called Hessian method, which has been discussed in detail 
in [143, 144]. Being the minimum X^{p) characterized by the set of parameter values p^, this 
method is based on the approximation of the variations of around its minimum by: 

AX' = x'M-xHp") ^ E ^A^, + -_£ (2.5) 

where d is the number of parameters considered in the error estimation (see below), Api=pi — p'- 
are displacements from their values at the minimum and the derivatives are also taken at the mini- 
mum, where the linear term in the above equation should vanish (meaning that the calculations of 
physical observables in should vary linearly around this point in the parameter space). There- 
fore all the information of x^ around its minimum is contained in the matrix of second derivatives 
or Hessian matrix Hjj^ in this approximation. The idea behind the method is to use this matrix to 
calculate the variations of the results with the parameters in the "vicinity" of the minimum, where 
the approximation should work. What is meant by "vicinity" is characterized by the tolerance 
parameter T via AX'^ < T^. 

If the data used were perfectly uncorrelated and perfectly compatible between them, the toler- 
ance parameter corresponding to 1 (7 uncertainty in the calculations (stemming, of course, from 
the 1 a experimental errors in x^) should be r= 1. As in others QCD fits (cf. [86, 143]), this is 
not our case because our data come from different observables and experiments and, moreover, 
the correlated errors are included in the experimental error in Eq. 2.4. As, for example, which 

* We define the Hessian matrix as Hij = , with no factors of one half. 
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deviations from the expected value of it y^lN are acceptable for 'X^ in a global QCD fit, which 
value to choose for the tolerance parameter T is somewhat arbitrary and subject to subjective 
interpretations. Since it seems reasonable that it somehow scales with the number of data (sets 



and) points included in the fits, we have chosen = (j^yi ; this results in T =4.654 (4.535) for 



our LO/NLO (NNLO) fits. Here we interpret the (1 cr) deviation (V^N) in as more than one 
standard deviation in the predictions, say 90%, and we have rescaled conveniently to what we 
consider as one standard deviation errors; again, this is completely arbitrary. 

For the numerical evaluation of the Hessian matrix it is convenient to map the parameter 
space into a coordinate system where the tolerated "vicinity" of the minimum is represented by a 
hypersphere. In this way the variations of X'^ are uniform and the instabilities resulting from the 
numerical difference of the typical changes in X'^ in different directions of the parameter space are 
avoided [144]. Since the Hessian matrix is symmetric, it has a complete set of (d) orthonormal 
eigenvectors V^- = vnc with eigenvalues e^, that is: 

d d d 

Y^HijVjk = ekVik, Y^VijVik = 5jk, then Hij = EkVikVjk (2.6) 

7=1 !=1 k=l 

The rescaled eigenvectors of the hessian matrix provide a natural basis for the mentioned trans- 
formation, which is achieved via: 

^ [2 
^Pi = E ^i^^^^ = \ ir'^ik, (2.7) 

k^i V ek 

where Zk are the new coordinates in term of which the tolerated region is the hypersphere: 

d 

AX^ = Y,zl<T\ (2.8) 

k=l 

In practice the matrix Mij is calculated by an extension of MINUIT called ITERATE [144] which 
uses an iterative procedure which converge to the eigenvectors and eigenvalues of the Hessian 
matrix of x^ in Eq. 2.4. 

With help of Mjj it is possible to calculate the eigenvector basis sets consisting in 2d sets of 
parameters pf-' = p° ± TMij which contain (in conjunction with T) the same information than 
the hessian matrix. Although by this construction the number of parameters entering in the error 
calculations is (unnecessarily) increased, it turns out to be useful for the propagation of experi- 
mental errors into quantities which depend on parton distributions. For example, the derivatives 
and a imit vector in the direction of the gradient of any quantity X(p) are: 

dX X{p+^)-X{p-^) X{p^^)-X{p-^) 

^ = 1^ , Dk[X) = — (2.9) 

2r ^i:j{x{p+i)-x{p-jY 

and we can estimate the uncertainty which follows from the condition /^x^ < by a displacement 
of length T in the direction of the gradient; it is simply given by: 
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Figure 2.2: Eigenvalues of the Hessian matrix for our LO, NLO (MS) and NNLO (MS) fits ordered from largest 
to smallest. In this order their decrease is, as expected [143], approximately exponential. 



Considering that the calculations required for practical applications in collider phenomenology 
are usually rather involved, the method is indeed quite convenient. Tables with the parameters of 
these eigenvector sets corresponding to our different fits are given in App. A; the eigenvalues of 
the hessian matrix (cf. Fig. 2.1) are also included. 

For the parameters entering in the error estimations the expressions above reduce to: 




(2.11) 



Since M/j-oc a large (small) eigenvalue corresponds to a direction in which increases 
rapidly (slowly), making the parameters tightly (poorly) constrained; they are usually referred 
to as steep (flat) directions. Note as well that the elements of the unit-length gradient indicate 
the influence of a particular eigenvector over the error of a quantity; in the case of the fitting 
parameters they give, by comparison, also information on their correlations. As suggested in 
[95], we include in our final error analysis only those parameters that are actually sensitive to the 
data used, i.e., those parameters which are not close to "fiat" directions in the overall parameter 
space. For the used data and our functional form Eq. 2.2 <i = 13 such parameters, including 
as{M^), are identified and are included in our final error analysis; the remaining ill-determined 
eight polynomial parameters Af and Bf, with highly-correlated uncertainties of more than 50%, 
were held fixed. 
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At this point it should be mentioned that we have explicitly checked that the linear approxi- 
mation is appropriate for our results, i.e., that the x'' variations around the minimum are approxi- 
mately quadratic for all parameters included in the error analysis. It is also worth to mention that 
since the value of T is rather arbitrary, in order to compare uncertainty results for distributions 
with different conventions one has to take into account the different definitions of the tolerance 
parameter. Fortunately, as is clear from its expression, the results for different conventions are 
linearly scalable, i.e., t^Ti= y.^Tj- 



2.4 The New Generation of Dynamical Parton Distributions 

In this section we compare our new dynamical parton distributions with the results obtained in 
recent common "standard" analyses, including the results of our "standard" fits which, as men- 
tioned, were carried out in exactly the same conditions as their dynamical counterparts. Recall 
also that in our analyses, as in all presently available distributions, the renormalization and factor- 
ization scales were set equal, which simplifies considerably the results but (obviously) prohibits 
their independent variation; we denote them here and in the following by = jX^ = jXp = jl^. 

Our dynamical input distributions are presented in Fig. 2.2 according to the parameters in 
Table. A.l. Being the dominant parameter in the small-x region Uf, the dynamical distributions 
have by construction, i.e., by optimally choosing a low /Iq, a valence-like structure (a/>0) for 




Figure 2.3: The valence-like input densities at /iy = 0.3, 0.5, 0.55 GeV^ for our dynamical LO, NLO (MS) and 
NNLO (MS) results respectively. The strange sea .s = 5^ vanishes at the input scale. Note that the lines of the gluon 
input distributions have been plotted thicker than the other lines in order to avoid confusion. The strange sea s = s 
vanish at the input scale. 
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all partons; in other words, not only the valence but also the sea and gluon input distributions 
vanish at small-^:. This is in contrast to the gluon and sea distributions of the "standard" fits, 
where the input scale is /Iq =2 GeV^ and thus have fl/<0 as is usual for common "standard" fits. 
Note that in both cases the valence distributions as well as A have a strong valence-like behavior, 
being in the former case enforced by the quark number constrains in Eq. 2.3 while in the latter it 
indicates that the small-x sea turns out to be approximately isospin symmetric. It should also be 
emphasized that all of our valence-like input distributions as well as the ones for the "standard" 
fits are manifestly positive in contrast to the negative gluon distributions in the small-^; region 
believed to be needed in other standard fits [84, 85, 107, 108, 110]. 

Our [3] NLO (MS) valence-like input distributions together with their 1 o uncertainties are 
compared with the previous GRV98 [1 19] NLO (MS) dynamical distributions at their respective 
input scale in Fig. 2.3. As can be seen the results of both analyses turn out to be very similar except 
for the gluon which peaks at a slightly larger value of x. However, such differences are merely 
within a 2a band of the new dynamical results. Both the LO and NLO (DIS) results are also 
similar in both analyses. Note however that the DIS distributions in [119] were not determined 
through a fit in the DIS factorization scheme but rather obtained through the factorization scheme 
transformations implied by Eq. 1.53. This treatment leads sometimes to a small negative gluon 
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distribution in the large x > 0.5 region which, however, is compatible within uncertainties with 
our definite positive gluon distribution resulting from a full analysis using the DIS factorization 
scheme. Besides this observation, the DIS results are, as mentioned, quite similar in both analysis. 
Furthermore by transforming our distributions from one scheme to the other we get also very 
similar results, which confirms the applicability of the transformations of Sec. 1.8. 

With respect to the uncertainties obtained in our different analyses, it must be noted that the 
dynamical input implies, as expected, a stronger constrained gluon distribution at larger values 
of as compared to the gluon density obtained from "standard" fits with a conventional non- 
valence-like input at /Iq > 1 GeV^, as can be seen in Fig. 2.4 for our NLO (MS) results. Since our 
valence-like sea input has a rather small a^, i.e., vanishes only slowly as ;c — > 0, the uncertainties 
of the sea distributions {ii + d) turn out to be only marginally smaller than those of the "standard" 
fit where the sea increases as x — > already at the input scale. Notice that the uncertainties 
generally decrease as increases due to the QCD evolution, as observed also in other analyses 
[86, 143]. 




Figure 2.5: The smaU-x NLO (MS) predictions of our dynamicaUy (radiatively) generated gluon and sea-quark 
distributions, together with their ±1 C7 uncertainties, as compared to the resuks of a "standard" analysis carried 
out under the same conditions. To ease the visibility of the two error bands of x{u + d) we have multiplied the 
stdNLO results by 0.9 as indicated. The corresponding GRV98 predictions lie within the 1 o band of our new 
dynNLO results. 
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Figure 2.6: Comparing the ± 1 cr error bands of our dynamical (dyn), standard (std) and CTEQ6 [95] NLO (MS) 
gluon distribution at small x. Note that Q^=2 GeV^ is the input scale of our "standard" fit, which is close to the 
CTEQ6 input scale ^Q=m^ ~ 1.7 GeV^, where the standard CTEQ6 fit employs a valence-like gluon input. The 
results at 2^ = 2 and 20 GeV^ have been multiplied by 0.5 and 2 respectively as indicated in the figure. 
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Figure 2.7: As Fig. 2.5 but for the sea quark distributions u + d. 
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Figure 2.8: Comparing our dynamical NLO (MS) ;/ and d parton distributions with the ones of CTEQ6.5 at 
2^ = 10 GeV^. The shaded areas represent the estimated ±1ct uncertainty of our analysis. 



At this point it should be mentioned that the standard CTEQ6 analysis [95] surprisingly re- 
sulted in a valence-like input gluon distribution at a scale as large as /Iq ~ 1.7 GeV^. Thus 
the CTEQ6 gluon distribution is expected to be similarly tightly constrained at (£■ > /Xq as our 
dynamical results. That this is indeed the case is illustrated in Fig. 2.5 where the 1 a uncertain- 
ties^ of the CTEQ6 gluon are similar in size to our dynamical NLO (MS) result, whereas the 
"standard" fit, with Q^ = 2 GeV^ and a non-valence-like structure, results in a sizeably larger 
uncertainty. Due to the sizeably different input scales the CTEQ6 gluon falls up to 30 - 40% 
below our dynamical NLO gluon for x< 10^^ and 2^ > 10 GeV^. 

The situation is different for the sea distribution in the small-x region (cf. Fig. 2.7); here the 
CTEQ6 input at /^q = ml is not valence-like, i.e. increases as x^O as expected for a "standard" 
fit, and thus the 1 cr uncertainty is comparable to our "standard" NLO fit, both being larger than 
the 1 a uncertainty obtained from our dynamical NLO fit. 

It should be emphasized that the distributions resulting from (GM)VFNS analyses (i.e., the 
newer sets of CTEQ and MRST) depend strongly on the model assumed for the treatment of 
heavy quark contributions; which is in contrast to FENS fits where heavy quark mass effects are 
fully taken into account in a theoretically unambiguous framework. For example, the inclusion 
of finite charm mass effects in CTEQ6.5 [98] reduces the charm contribution to F2{x, Q^), which 
is compensated by larger u and d distributions at small x as compared to the standard CTEQ6 
analysis [95], where heavy quarks have been treated in the zero-mass approximation. Fig. 2.7 
illustrates that such "enhanced" quark distributions differ very little from our dynamical results. 
These ratios remain practically unchanged at higher scales, like =M^ relevant for pro- 
duction, therefore our predictions for hadronic W^/Z^ production cross-sections, for example, 
at Tevatron and LHC are similar to the "enhanced" ones observed in [98]. We study hadronic 
weak gauge boson production in more detail in Sec. 2.10. 



'when comparing our uncertainty results with the ones of CTEQ6 [95] where T = 10 has been assumed, we rescale these CTEQ 
uncertainties according to A«, 0.4654Afli in order to comply with our T = 4.654. 
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Figure 2.9: Comparing the present dynamical NLO (MS) distributions with the previous dynamical ones of 
GRV98 at = 10 GeV^. The shaded areas represent the estimated ±1 CT uncertainty band of our analysis. 

Moreover, this holds also true for the GRV98 analysis [119], indeed in the relevant small-x 
region these ratios would be practically unaltered if the GRV98 distributions were used instead of 
our current result since both dynamical analyses give rather similar results. Despite the fact that 
their respective analyses were done with approximately ten years of time difference and mean- 
while a substantial increase in the precision of the HERA data has been achieved, the dynamical 
distributions, i.e. the predictions stemming from the dynamical model, remain stable to within 
less than about 20%. This is illustrated in Fig.2.8 where ratios of our dynamical NLO (MS) gluon, 
sea and valence distributions are compared with the previous GRV98 results and at a represen- 
tative scale Q^ = 10 GeV^ relevant for phenomenological applications. As we have seen, this is 
in contrast to "standard" fits, which depend strongly on the particular details of the analyses (e.g. 
the data included, ansatz used, details on the treatment of quark masses, etc.). 

Turning now our attention to our [5] NNLO results, we present in Fig. 2.9 our dynamical 
NNLO (MS) input distributions at = 0.55 GeV^, according to the parameters in Table A. 1 , to- 
gether with their 1 a uncertainties; our NLO (MS) input distributions are shown again for compar- 
ison. The resulting valence {Uy and dy) distributions are now somewhat enhanced around x = 0.l 
to 0.2 respect to our NLO (MS) results, and a strong and clear valence-like small-^: behavior of 
the NNLO gluon input is now required, ag = 0.994ib0.379, as compared to Og ~ 0.5168± 0.4017 
at NLO. Furthermore there is also a strong enhancement of the NNLO gluon over the NLO one 
around x = 0.1 and a sizeable depletion at larger values of x. The 1 o uncertainties at NNLO are 
comparable to the NLO ones shown in Fig. 2.3, except for the NLO gluon at x > 0.3 which is 
stronger constrained due to the light high-pT jet data (cf. Sec 2.2) which cannot be consistently 
included in an NNLO analysis. 
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Figure 2.10: The dynamical NNLO(MS) input distributions and their ±1 a uncertainty band at /Xq =0.55 GeV^ 
together with our dynamical NLO (MS) input distributions at /Iq = 0.5 GeV^. The strange sea s = s vanishes at 
the input scale. 

The valence distributions of our NNLO "standard" analysis (jUq = 2 GeV^) are compared with 
the "standard" NNLO ones of Alekhin, Melnikov and Petriello [105] (AMP06, Mo =9 GeV^) 
and the "standard" NNLO pure NS analysis of Blumlein, Bottcher and Guffanti [147] (BBG06, 
/^o = 4 GeV^) in Fig. 2. 10 at 2^ = 4 GeV^. In the relevant valence x-region, x > 0. 1 , we confirm 
the BBG06 results, in particular the enhancement of x<i,, with respect to the result of AMP06. In 
any case we, as well as BBG06, observe a significant enhancement of the NNLO xuy and xdy 
with respect to the NLO results. Our dynamical valence distributions at =4 GeV^ practically 
coincide with the standard ones shown in Fig. 2. 10, which indicates that the valence distributions 
are very robust with respect to the choice of the input scale /Iq. 

Again, the distinctive valence-like NNLO gluon input at low I^q <l GeV^ of the dynamical 
model (cf. Fig. 2.9) implies a far stronger constrained gluon distribution at larger values of 
as compared to a gluon density obtained from a NNLO "standard" fit with a conventional non- 
valence-like input at /^Q > 1 GeV^, namely I^q=2 GeV^, as can be seen in Fig. 2.1 1. 

In contrast to our "standard" NLO results, the "standard" NNLO gluon input at jUq =2 GeV^ 
is very weakly constrained at small x (flg = 0.0637 ±0.1333, cf. Table A.l) and therefore mo- 
mentum conservation (cf. Eq. 2.3) cannot sufficiently constrain it at larger values of x, moreover 
since high-z^r jet data have not been taken into account as already pointed out. Note as well that 
this common "standard" NNLO input gluon distribution at /Xq =2 GeV^ is also compatible with a 
valence-like small-x behavior {ag>0), a tendency already observed in [85], and that our dynam- 
ical NNLO gluon distribution in Fig. 2.1 1 remains valence-like, i.e. decreases with decreasing x. 
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Figure 2.11: Our "standard" NNLO(MS) valence distributions and their ±lcT uncertainties at Q^=4 GeV^ 
together with our "standard" NLO(MS) results and the "standard" NNLO results of AMP06 [105] and BBG06 
[147]. Our dynamical valence distributions at 2^ = 4 GeV^ practically coincide with the standard ones shown. 

even at Q^ = 2 GeV^. This is mainly caused by the NNLO splitting function P^g^ (cf. Sec. 1.2 to 
1.4) which is negative and more singular in the small-.x region [196] than the LO and NLO ones^. 

As in the NLO case, the uncertainties generally decrease as increases due to the QCD 
evolution, but the ones of the dynamical predictions in the small-;*: region remain substantially 
smaller than the uncertainties of the common "standard" results as exemplified in Fig. 2.11. Fur- 
thermore, it is a general feature of any NNLO gluon distribution in the small-;c region that it falls 
below the NLO one as can be seen in Fig. 2.11 by comparing the solid curves with the long- 
dashed ones, and the dashed-dotted curves with the dotted ones. In the first dynamical case the 
NNLO predictions for x < 10^^ are several (7 below the NLO ones. 

The dynamical sea distribution x{U + d) derives from a less pronounced (az < ag) valence-like 
input in Fig. 2.9 which vanishes very slowly as a:— >0; <3j; = 0.1374ib0.0501 (cf. Table A.l). This 
implies that the valence-like sea input is similarly increasing with decreasing x down to ;c ~ 0.01 
as the sea input obtained by the common "standard" fit where = —0.1098 ±0.0122 according 
to Table A. 1 . Therefore, the 1 o uncertainty bands of our dynamically predicted sea distributions 
at larger values of in Fig. 2.13 are only marginally smaller than the corresponding ones of 
the "standard" fit. In contrast to the evolution of the gluon distribution in Fig. 2.11, the NNLO 
sea distributions in Fig. 2.13 lie always above the NLO ones in the small-x region, x< 10^^, and 
at not too large values of Q^. Here all NNLO sea distributions are rather similar, including the 
"'standard" one of A1VIP06 [105] (cf. Fig. 2.13). 

*For fi/ = 3, -vPjg' (.v) -3147.66 In | + 14737.89 as x 0, whereas xP^f {x) ^ 12 and jcpj]' (x) ^ -81 .33. 



42 



2.4. The New Generation of Dynamical Parton Distributions 




Figure 2.12: Comparing the ±lcT error bands of our dynamical (dyn) and standard (std) NNLO(MS) gluon 
distributions at small x for various fixed values of . Our NLO (MS) dynamical and "standard" results as well 
as the "standard" NNLO results of AMP06 [105], based on an input scale ^q = 9 GeV^, are shown as well. Note 
that Q^ = 2 GeV^ is the input scale of our standard fits. The results at 2^ =2 and 20 GeV^ have been multiplied 
by 0.5 and 2 respectively as indicated in the figure. 
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Figure 2.13: As in Fig. 3 but for the sea quark distribution x{u + d). 
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2.5 Determination of the Strong Coupling 

As already mentioned, we use the reference value as{M^) as contour condition for the determina- 
tion of the strong coupUng constant; for which we always employ the exact iterative solutions (cf . 
Sec. 1.1). Instead of using a fixed value for as{M^), typically the LEP measurement or the world 
average of 0.1176 [19], we include it as a free parameter in our fits and their error estimations. 
Thus the running coupling is determined in our analyses together with the parton distributions 
of the nucleon, in particular it is closely related to the gluon distribution which drives the QCD 
evolution. We briefly discuss in this section the results obtained and compare them with similar 
determinations of the strong coupling. For a more detailed and comparative recent discussion of 
NLO and NNLO results the reader is referred to [148]. 

It is worth to mention at this point that our treatment implies that each eigenvector basis set 
(cf. Sec. 2.3) has a different associated value of as{M^)- The use of these values is mandatory for 
uncertainty studies; note that the difference between values of as{Q^) corresponding to different 
eigenvector sets (belonging to the same fit) can be as large as 10% at low scales. 

In Fig. 2. 13, we show the results for the running coupling as{Q^) according to the values 
obtained for as{M^) in our dynamical fits (cf. Table A. 1) together with their respective ±1 cr error 
bands. Note, as is well known, that the LO solutions are much larger and the results stabilize at 
NLO. In general, the higher the perturbative order the faster (Xs{Q^) increases as decreases. In 
order to compensate for this increase, a NNLO analysis is expected to result in a smaller value for 
as(M|) than a NLO one as is indeed the case. Both values lie, however, within a 1 a uncertainty. 
Similar results were obtained in a previous dynamical fit [149] which, however, was performed 
for a restricted set of (mainly small-x) DIS data. 

0.55 p-n-r| 1 I 1 I 1 I 1 I ' I ' : 




0.05 ^ — — — — — — ^ 

1 10 100 10^ 10"^ 10^ 



Figure 2.14: Values of (Xs{Q^) resulting from our dynamical fits together with their respective ±1 CT error bands. 
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According to the 1 a errors for as{M^) quoted in Table A. 1, in general the strong coupling 
is stronger constrained in the dynamical fits than in their "standard" counterparts. This is in 
principle expected due to the larger evolution distance of the dynamical distributions. It should 
also be mentioned again that there is a certain correlation between the chosen value of and the 
resulting values for aj(M|), which increase for increasing i^q. 

At NNLO the a^-uncertainty of our dynamical fit is about half as large as the "standard" 
one. Apart from the larger evolution distance, this is due to the strongly constrained valence- 
like input gluon distribution in the small-x region; consequently the energy-momentum sum rule 
(cf. Eq. 2.3) sufficiently constrains xg{x,[XQ) in the medium to large x-region in the dynamical 
case while in the "standard" one the gluon distribution, and thus a,(M|), is considerably less 
restricted, moreover since jet data were not included at NNLO. 

Keeping in mind that our stated errors always refer to 1 (7 uncertainties, our "standard" NLO 
fit error for o;v(M|) is compatible with the 2a uncertainty stated in the literature (see, e.g., [95] and 
the discussion in [84]). Our NLO (MS) results are also compatible, within about 1 <J uncertainty, 
with the ones obtained from analyzing only DIS structure functions. The same holds true for 
our "standard" NNLO fit result (cf. Table A.l) which, within errors, is compatible with the one 
obtained from a "standard" fit [150] to a restricted set of small-x DIS data. 

Previous "standard" NNLO fits considering only the flavor non-singlet valence sector of struc- 
ture functions [147, 151, 152] resulted in somewhat smaller values of as(M|) but remain within 
1 to 2a uncertainty. A similar NS valence analysis [153] as well as a full analysis [154] being 
based, however, on incomplete calculations of the moments of 3-loop anomalous dimensions 
(splitting functions) yielded slightly larger values of at NNLO, however the estimated errors 
are large enough so as to comply with our result in Table A. 1 . 

As already mentioned several times, data for high-pr jet production in hadron-hadron scatter- 
ing should not be included in a consistent NNLO analysis, since such processes are theoretically 
known only up to NLO (Sec. 1.10). Including them nevertheless in a "standard" NNLO analysis 
requires [88, 89] generally larger values of as{M^). 

Our ai(M|)-uncertainty in Table A.l is also about twice as large as the one obtained in a 
comparable standard NNLO analysis [105] where the high-pr jet data have been disregarded 
for consistency reasons as well. Without these data the gluon distribution is little constrained in 
the medium- to large-x region where it plays an important role for the QCD evolution at small 

(k) 

values of x due to the convolution with the dominant Pgg . This a^-uncertainty remains sizeable 
irrespective of the choice of the input scale Mo > ^ GeV"^. Only within a Bayesian treatment of 
systematic errors, by taking into account point-to-point correlations [155-157], the uncertainty 
of aj(M|) turns out to be about two times smaller [103-105, 155-157]. 

2.6 Comparisons with the Data Used 

As discussed in Sec. 2.2, we have used for our analyses the experimentally directly measured "re- 
duced" DIS cross-section a^^^ (cf. Eq. 1.41), although fitting just to the dominant (one-photon 
exchange only) structure function fJ resulted in very similar results. Representative comparisons 
of our results with the relevant HERA (HI and ZEUS) data [106, 123-126] on this structure func- 
tion of the proton, /^/(x, Q^), are presented in Figs. 2. 14 (NLO) and 2. 15 (NNLO). As can be seen 
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Figure 2.15: Comparison of our dynamical (dyn) and standard (std) NLO (MS) results for F2{x, Q^) with HERA 
data for > 1.5 GeV^ [106, 123-126]. The dynamical LO results are shown for comparison. To ease the 
graphical presentation we have plotted F^(x,2^) +0.5 x i{Q^) with i{Q^) indicated in parentheses in the figure. 
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2. Analysis and Results 



the data are well described throughout the whole medium- to small-;c region for 2^ > 2 GeV^ 
and thus perturbative QCD is here fully operative. It should be reemphasized that, due to our 
valence-like input, the dynamical small-;c results {x< 10^^) are predictions being mainly gener- 
ated by the QCD evolution; this is in contrast to a common "standard" fit where the gluon and 
sea input distributions do not vanish as x — > 0. At <2 GeV^ the theoretical results fall below 
the data in the very small-x region, although the NNLO results are closer to the data than the 
NLO ones; this is not unexpected for perturbative leading-twist results, since nonperturbative 
(higher-twist) contributions to F2 will eventually become relevant, even dominant, for decreasing 
values of Q^. 

Furthermore, the more restrictive ansatz of the valence-like input distributions at small-x as 
well as the sizeably larger evolution distance (starting at /^q < 1 GeV^) generally imply, as for the 
distributions themselves (cf. Sec. 2.4), smaller uncertainties concerning the behavior of structure 
functions in the small-x region than the corresponding results obtained from the common "stan- 
dard" fits, in particular as increases. This is illustrated in Fig. 2.16 for our NLO (MS) results 
where the error bands correspond to a 1 (7 uncertainty. However, since our NLO valence-like sea 
input has a rather small power of x, i.e. vanishes only slowly as x— >0 as already discussed, the 
uncertainties of the sea dominated F2{x, Q^) turn out to be not too different from the "standard" 
fit where the sea increases as x— >0, negative power of x, already at the input scale IJ.q=2 GeV^. 
The results for our dynamical and standard NNLO fits are respectively rather similar to the NLO 
ones shown in Fig. 2.16 (cf. Fig. 2.17). 




10"^ lO""^ 10"^ 0.01 0.1 1 

X 

Figure 2.17: Typical ±1 CT uncertainty bands of our dynamical and "standard" NLO (MS) results for two repre- 
sentative values of Q^. To ease the visibility of the two error bands at = 2.5 GeV^ we have subtracted 0.2 from 
the stdNLO result as indicated. For illustration two HI and ZEUS data points from Fig. 2.14 with their almost 
invisible errors are shown as well at 2^ = 2.5 and 15 GeV^ respectively. 
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Although the inclusion of NNLO corrections implies an improved value of j}, as mentioned 
typically Znnlo — ^-^/t^NLO ('^f- Table A.l), present high precision DIS data are not sufficiently 
accurate to distinguish between the NLO results and the minute NNLO effects of a few percent. 
This is illustrated in Fig. 2.17 where the experimental (statistical and systematic) errors are far 
bigger than the differences between the NLO and NNLO results. It should, however, be noticed 
that the NNLO uncertainty bands are somewhat narrower (reduced) than the ones at NLO; the 
results are similar for our "standard" fits. By analyzing only the flavor non-singlet valence sector 
of structure functions, it was already been pointed out [152] that NNLO effects cannot be delin- 
eated by present data in the medium- to large-^: region, and moreover, uncertainties of NLO and 
LO analyses (such as higher twists, different factorization schemes and QED contributions to the 
QCD evolution) turn out to be comparable in size to the NNLO 3-loop contributions [152]. 

Turning now to heavy quark production in DIS, recent HERA data on charm [134—137] and 
bottom [136, 137] contributions to f2(^, 2^) are compared in Fig. 2.18 with our dynamical LO 
and NLO (MS) predictions. Although not shown, it is worth to mention that also in this case the 
uncertainties of the "standard" results turn out to be larger than the ones implied by the dynamical 
model. The impressive agreement observed in Fig. 2.18 illustrates the reliability of the the nf = 3 
FFNS, in which the heavy quarks (h = c,b,t) are always produced as final states in fixed-order 
perturbation theory via hard production processes initiated by the light partons of the nucleon 
{u,d,s quarks and the gluon g). In other words, the intrinsic heavy quark content of the nucleon 
(if any) is marginal since all the observed heavy quark contributions are well taken into account by 
extrinsic heavy quark production, i.e., originated only from light quark partons in the initial state. 
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Figure 2.19: The dynamical NLO (MS) predictions for ' in the strict nf = 3 FFNS, choosing )l = 4m| with 
mc = 1.3 GeV and =4.2 GeV, together with the ±1 (T uncertainty band. For comparison we also display the 
central LO predictions which are entirely due to the 7*-gluon fusion subprocess y*g^ hh. The charm production 
data as obtained from D* measurements are taken from [134, 135] (solid and open squares) and the HI direct 
track measurements from [136, 137] (open circles); the bottom production data taken from [136, 137]. 



Therefore only the nf = 3 light u,d,s quark flavors and gluons constitute the genuine (intrinsic) 
partons of the proton and (in principle) the heavy c, b, t quark flavors should not be included in the 
parton structure of the nucleon. This fact together with the perturbative stability of heavy quark 
production (cf. Sec. 1.9), even at observed also in Fig. 2.18, indicate that there is no 

need to resum supposedly "large logarithms" (in^j); which is of course in contrast to genuine 
coUinear logarithms appearing in light (massless) quark and gluon hard scattering processes. In 
our opinion the use of the VFNS should therefore be avoided when possible; somewhat dissenting 
views are summarized in [158]. We will have a closer look at the above mentioned perturbative 
stability of heavy quark production in DIS and determine under which conditions the use of the 
(Z]VI)VFNS is appropriate in Sees. 2.9 and 2.10. 

As already pointed out (cf. Sec. 1.10), the measurements of Drell-Yan dilepton production 
in pp and pd collisions [138, 139] are instrumental in fixing A = d — u (or d/u). In Fig. 2.19 
we display our dynamical NLO and NNLO results, together with the ztlcr uncertainties, and 
the NNLO A1VIP06 ones from [105] for the differential dimuon mass distributions for pp and 
pd collisions respectively; the "standard" fit results differ only marginally from the dynamical 
ones. The data are given for various average values of xp =x\ — X2, where x\ and X2 refer to 
the fractional momenta of the quarks in the beam {p) and the nucleon target (A^) respectively; 
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Figure 2.20: Our dynamical NNLO and NLO [3] results, together with their ±1 CT uncertainties, for Drell-Yan 
dilepton production in pp and pd collisions respectively for various selected average values of Xf using the 
data sets of [138]. For comparison the NNLO AMP06 results [105] are shown as well. To ease the graphical 
presentation we have multiplied the results for the cross-sections by 10' with / indicated in parentheses in the 
figure for each fixed average value of Xf . 
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experimentally xf>0 {xi >X2) and consequently the Drell-Yan cross-section is dominated by 
the annihilation of a beam quark with a target antiquark. In the relevant kinematic region where 
high-statistics data exist, all three NNLO and NLO results shown agree within 1 a. 

In order to test the dependence of our results on the specific choice of the factorization scheme, 
we have repeated our dynamical NLO analysis in the DIS factorization scheme of Sec. 1.8. Since 
the (light) parton distributions in the DIS scheme are defined via the F2 structure function, it is not 
very surprising that, in contrast to the DIS parton distributions themselves, the results for physical 
observables directly related to DIS structure functions are very similar. Indeed the results in the 
DIS factorization scheme for F2{x,Q^) in Figs. 2.14 and 2.16 are practically indistinguishable 
from the ones in the MS scheme, as are the results for F2^{x,Q^) in Fig. 2.18. Differences 
become visible only for processes which are not directly related to DIS structure functions, such 
as Drell-Yan cross-sections, but there the differences lie well within the ± 1 (7 uncertainty of 
the NLO (MS) results as illustrated in Fig. 2.20 for the DY-asymmetry. There we display our 
dynamical NLO (MS) result for o^'^ jloPP together with the it 1 a uncertainty band as well as the 
previous GRV98 result which agree in the statistically relevant A:-region, with x^ referring to the 
average fractional momentum of the target partons. The dynamical NNLO (MS) as well as the 
results from the "standard" NNLO AMP06 fit [105] are only marginally different (cf. [5]); being 
the differences between all the theoretical predictions much smaller than the experimental errors 




Figure 2.21: Our dynamical NLO result in the MS factorization scheme, together with its ±1 CT uncertainty, for 
with = ^j^j^f^, as a function of the average fractional momentum X2 of the target partons. The GRV98 
NLO result [119] is shown for comparison. The dynamical NLO (DIS) result is shown by the dashed-dotted 
curve. The data for the dimuon mass range 4.6 < < 12.9 GeV are taken from [139]. 



52 



2.7. Astrophysical Relevance of the Dynamical Predictions 




100 200 300 400 500 600 

Pt (GeV) 

Figure 2.22: The pp CDF high-px inclusive jet data [141] compared with our dynamical LO and NLO(MS) 
results, as well as with the NLO CTEQ6 result. 

and than the 1 a uncertainty bands obtained in our NLO and NNLO analyses; which differ only 
marginally as well (cf. [5]). 

Finally the pp CDF (Tevatron II) high-;?! inclusive jet data [141] are compared in Fig. 2.21 
with our dynamical LO and NLO (MS) results, as well as with the ones of CTEQ6 [95]. The 
small 1 (7 error bands are almost invisible on the huge logarithmic scale used. Our NLO result 
almost coincides with the one of CTEQ. There is a clear improvement at NLO as compared to LO 
which falls slightly below the data at 300 GeV. Nevertheless the LO YAgh-pi fit corresponds 
to ;|^^/dof ~ 1, which is only twice as large as at NLO (cf. Table A.l). 



2.7 Astrophysical Relevance of the Dynamical Predictions 

The predictive power of the dynamical model is especially important for investigations concern- 
ing vA^ cross-sections [160-162] of ultrahigh energy cosmic neutrinos produced via astrophysi- 
cal processes, e.g. in active galactic nuclei, black holes or in the decays of very massive particles 
(see, for example, [163, 164] for more details). Here one needs a somewhat reliable knowledge 
of parton distributions at the weak scale =M^ down to 10^^ (;c~ 2m^Ev ) highest 
energies of Ey ~ 10^^ GeV which in the "standard" approach requires extrapolations into yet un- 
measured regions. As our parameter-free small-^: predictions for parton distributions at x < lO^'^ 
are entirely of QCD-dynamical origin and depend rather little on the detailed input parameters 
at x> 10^^, it is interesting to study these predictions in kinematic regions not accessible by 
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Figure 2.23: Dynamical NLO (MS) predictions for F2{x, Q^), together witii their 1 CT uncertainties, for extremely 
small values of x. The 1 a uncertainties of the standard (std) NLO fit extrapolations are shown by the vertical 
bars. The dotted curves are the contributions from the light (np = 3) quarks and gluons to Fj' for the dynamical 
(dyn) NLO result. In other words, the difference between the dotted and solid curves is due to NLO heavy quark 
(charm, bottom) contributions which derive from photon-gluon (quark) fusion processes. The dynamical GRV98 
predictions lie within the ±1 CT band of our present dynNLO predictions. 

present DIS experiments. Furthermore, the more restrictive ansatz of the valence-like input dis- 
tributions at small-.x as well as the sizeably larger evolution distance imply, as we have seen, 
smaller uncertainties concerning the behavior of structure functions in the small-x region than 
the corresponding results obtained from the common "standard" approach, in particular as 
increases. 

Since Fj^x^Q^) is, in the very small-jc region, dominated by u{x,Q^) and d{x,Q^), as are 
the CC neutrino-(isoscalar) nucleon cross-sections, the F^' structure function can be utilized 
for estimating the magnitude of uncertainties of the predictions in the extreme small-.x region, 
which are shown in Fig. 2.22. At =^w dynamical NLO predictions correspond to a ±1 cr 
uncertainty of about ±7% at;c= 10^^ whereas the uncertainty of the extrapolation of a "standard" 
fit is about twice as large. At smaller scales the uncertainties obviously increase as illustrated in 
Fig. 2.22 at = 100 GeV^. Taking into account previous extrapolation ambiguities [119], one 
can conclude [160] that the dynamically predicted small-.x parton distributions allow neutrino- 
nucleon cross-sections to be calculable with an accuracy of about 10% at the highest cosmic 
neutrino energies. It should be mentioned that an ad hoc fixed power law of x extrapolation of 
the standard CTEQ6.5 structure functions [98] to x= 10^^ at Q^=M^ lies, accidentally, only 
about 10% below our dynamical NLO prediction in Fig. 2.22. On the other hand, an alternative 
parametrization [165] of present HERA (ZEUS) data which is not QCD oriented but based on 
analyticity and unitarity gives, when extrapolated to x= 10^^, a factor of about 6 smaller a value 
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forf/(10"^M^) than our result in Fig. 2.22. Since the perturbative dynamical QCD /predictions 
for the small-A; behavior of structure functions down to x= 10^^ proved to be in agreement with 
later HERA measurements as discussed in Sec. 2.6, it is hard to imagine that perturbative QCD 
dynamics and evolutions should become entirely inappropriate atx= 10^^ to 10^^ at even much 
larger scales. 



2.8 The Longitudinal Structure Function 

The importance of using the experimentally directly measured "reduced" DIS cross-sections a^'^ 
(cf. Eq. 1.41) for the determination of parton distributions has been emphasized in [110]. As 
mentioned, for not too large values of this quantity is dominated by the one-photon exchange 
structure functions, in particular by F2. However, the effect of Fl becomes increasingly relevant as 
X decreases at a given Q^, where y increases. This is seen in the data as a flattening of the growth 
of Or{x,Q^) as X decreases to very small values, at fixed Q^, leading eventually to a turnover as 
can be seen in Fig. 2.23. At the lower values of in Fig. 2.23 it was not possible in [110] to 
reproduce this turnover at NLO. This was mainly due to the negative longitudinal cross-section 
(negative Fl{x,Q^)) encountered in [85, 110]. Since all of our cross-sections and subsequently 
all structure functions are manifestly positive throughout the whole kinematic region considered. 
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Figure 2.24: The dynamical NNLO predictions, together with their ±1 a uncertainties, for the "reduced" DIS 
cross-section Orix, Q^) = F2 — (y^ /Yj^)Fi. The uncertainty bands of our previous dynamical NLO results (dashed 
curves) are very similar in size [3] as the ones shown for NNLO. The HERA data for some representative fixed 
values of are taken from [106, 123-126]. 
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Figure 2.25: Our dynamical NNLO predictions for (Jr{x,Q^) for different proton beam energies Ep relevant 
for most recent HERA-Hl measurements [166]. The ±1 (7 uncertainty bands are similar to the ones shown in 
Fig. 2.23. Notice that the curves terminate when y = I. 

our dynamical NLO and NNLO results in Fig. 2.23 are in good agreement with all small-^: HERA 
measurements [106, 123-126]. The same holds true for our "standard" NLO and NNLO results 
which, besides having wider uncertainty bands, are almost indistinguishable from the dynamical 
ones shown in Fig. 2.23. In Fig. 2.24 we display our NNLO results for the "reduced" cross- 
section at different proton beam energies Ep, relevant for most recent HI measurements [159]. 
The turnover at small x becomes more pronounced at smaller energies because of the larger 
values of y. Our dynamical small-;c predictions are fully compatible with the (preliminary) HI 
data presented in [166]. 

Turning now to Fl itself we present our dynamical and "standard" LO, NLO and NNLO 
predictions Fig. 2.25, together with the small subdominant charm contributions at LO and NLO. 
These gluon-dominated dynamical predictions become perturbatively stable already at > 2-3 
GeV^, where precision measurements could even delineate NNLO effects in the very small- 
X region. This is in contrast to the common "standard" results, as has been already observed in 
[149], but the differences between the NNLO and NLO results are here less distinguishable due to 
the much larger 1 a uncertainty bands which partly overlap in the very small-^: region. It should 
be noticed that the NLO/NNLO instabilities implied by the "standard" fit results obtained in 
[85, 110] at e2<5GeV2 are far more violent than the ones shown in Fig. 2.25 for our "standard" 
results, which is mainly due to the negative Fl{x,Q^) encountered in [85, 110]. The perturbative 
stability in any scenario becomes in general better the larger Q^, typically beyond 4 — 5 GeV^ [36, 
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2.8. The Longitudinal Structure Function 





Figure 2.26: Dynamical and "standard" parton model LO, NLO and NNLO predictions for Fi{x,Q^) together 
with the ±1 CT uncertainty bands at NLO and NNLO. The heavy charm (c) contributions at LO (short-dashed 
curves) and NLO (long-dashed curves) are shown as well. The results at = 100 GeV^ are multiplied by 0.5 
as indicated. 



57 



2. Analysis and Results 



10"- 



10- 



10" 



10" 



1 

0.8 
_ 0.6 

a 

^ 0.4 
0.2 




— 1 1 1 1 1 — I — I — r- 



W = 276 GeV 



— 1 1 1 — I — I — r- 



— 1 1 1 1 1 — I — I — r- 



dynNNLO 

dynNLO 

stdNNLO 

stdNLO 



prelim. • 
e"'"p ■ 
e"p □ 




_i I I I I I I i_ 



_i I I I I I I i_ 



_i I I I I I I i_ 



10 



100 



10-^ 



q2 (GeV^) 



Figure 2.27: Dynamical and "standard" NNLO and NLO predictions for Fi{x, Q^) at a fixed value of W = 276 
GeV. The NLO (MS) results are taken from [3]. The (partly preliminary) HI data [106, 125, 167, 168] are at fixed 
W ^ 276 GeV. 



85, 110, 149], as evident from Fig. 2.25. This is due to the fact that the evolutions eventually 
force any parton distribution to become sufficiently steep in x. It should be mentioned that the 
sizeable discrepancies between NNLO and NLO predictions at Q^ = 2 GeV^ and xc^ 10^^ in 
Fig. 2.25 are not too surprising since 2^ ^2 GeV^ represents somehow a borderline value for the 
leading twist-2 contribution to become dominant at small-^: values. This is further corroborated 
by the observation that the dynamical NNLO and NLO twist-2 fits slightly undershoot the HERA 
data for F2 at <2 GeV^ in the small-^; region (cf. Figs. 2.14 and 2.15), which indicates that 
nonperturbative (higher-twist) contributions to structure functions become relevant there [3, 1 19]. 

For completeness we finally compare in Fig. 2.26 our (leading-twist) NNLO dynamical and 
"standard" predictions for Fl{x, Q^), together with their ±1 (7 error bands, with a representative 
selection of (partly preliminary) HI data [106, 125, 167-169] at fixed W ~ 276 GeV. For com- 
parison we also show in Fig. 2.26 our NLO results which have comparable, although somewhat 
larger, uncertainty bands (cf. [3]); note as well that the uncertainties of the dynamical predictions 
are again smaller than the corresponding to the "standard" fits. All our NNLO and NLO results 
for Ft, being gluon dominated in the small-.x region, are in full agreement with present measure- 
ments which is in contrast to expectations [84, 85, 110] based on negative parton distributions 
and structure functions at small values of x. To illustrate the manifest positive definiteness of our 
dynamically generated structure functions at >Mo' show Fi^{x,Q^) in Fig. 2.26 down to 
small values of although, as mentioned, leading twist-2 predictions need not necessarily be 
confronted with data below, say 2 GeV^. 



58 



2.9. The Role of Heavy Quark Flavors in DIS 

2.9 The Role of Heavy Quark Flavors in DIS 

In this and the following section we compare fixed flavor number scheme (FFNS) and (zero- 
mass) variable flavor number scheme (VFNS) parton model predictions at high energy colUders 
at LO and NLO of QCD. Based on our LO- and NLO-FFNS dynamical and "standard" parton 
distributions [3], we generate radiatively sets of VFNS parton distributions [4], where also the 
heavy quark flavors are considered as massless partons within the nucleon. In this section we 
focus on the role of these distributions in deep inelastic scattering, while in the next one we 
study their implications for hadron colliders. In particular we show that the VFNS predictions 
are compatible with the reliable FFNS ones for processes in which the invariant mass of the 
produced system far exceeds the mass of the participating heavy quark flavor [4]. 

As explained in Sec. 1 .9, the VFNS is characterized by increasing the number of massless 
partons np by one unit at the unphysical thresholds = m\, starting from = 3 at = '"c' 
i.e. c(x, m^) = c(;c, m^) = up to NLO. Thus the "heavy" n/r > 3 quark distributions are perturba- 
tively uniquely generated from the n^r — I ones via the massless renormaUzation group equations 
(evolution). The VFNS heavy quark contributions are therefore calculated in analogy to those of 
the light quarks, e.g. their contributions to DIS structure functions follow Sec. 1.7. Our choice 
for the input of the "heavy" h = c,b VFNS distributions are our FFNS distributions at = m\. 
The VFNS predictions at scales Q^^ml should become insensitive to this, somewhat arbitrary, 
input selection [103], whose virtues are the fulfillment of the standard sum-rule constraints (cf. 
Eq. 2.3) together with reasonable shapes and sizes of the various input distributions. 

Since the deep inelastic structure function F2{x, Q^) plays an instrumental role in determining 
the parton distributions of the nucleon, we consider as a first test of the VFNS "heavy" quark 
distributions the charm and bottom contributions to this quantity. In Fig. 2.27 we compare our 
VFNS results, with = as usual^ (cf. Sec. 1.7), with recent HERA data on heavy quark 
electroproduction. It can be seen that the data are, except for higher Q^, very badly described by 
the VFNS results. This well-known fact is moreover totally expected since the VFNS appear as 
a resummation of the FFNS contributions, which make sense only for S>m^. For this same 
reason it is not common to find comparisons of (ZM) VFNS predictions with heavy quark pro- 
duction data, especially near threshold; we show them here nevertheless because while presenting 
our results we got involved in discussions where we realized that the reader may not be aware 
of (or even reluctant to accept) this well-established fact. Note as well the excellent agreement 
of the FFNS predictions (already discussed in Sec. 2.6) and that, despite the common belief that 
"non-coUinear" logarithms (in ^) need to be resummed, for Q^^m\ present data are not suf- 
ficiently accurate to discriminate between VFNS and FFNS results in Fig. 2.27. At this point it is 
worth to recall that, in contrast to what it is sometimes assumed, the calculations in both schemes 
do not converge in the limit ^ since, even in this limit, the FFNS results always contain 
contributions from the (physical) threshold region 5th=4m| (cf. Sec. 1.9 and [48]). 

'Notice that here and in the following /i = jUp = /ij where Hr and ^lF are the renormaUzation and factorization scales respectively. 
As mentioned, this choice is dictated by the fact that our (and all other presently available) parton distributions were determined and 
evolved with ^f = Hr, i.e. with the commonly adopted standard evolution equations. 
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Q2(Gev2) Q^GeV^) 

Figure 2.28: Comparing our dynamical NLO(MS) calculations for {x,Q^) in the (zero-mass) VFTSTS with 
heavy quark electroproduction data [134-137]. The NLO(MS) FFNS predictions are shown for comparison; 
here /i^ = 4m|, with m^j, = 1.3,4.2 GeV, although the results are rather unaffected by this choice [48]. 

In Figs. 2.28 and 2.30 we compare the VFNS (/i^ = Q^) with the FFNS predictions for F:^{x, Q^) 
and f2^(x, Q^) respectively, here we choose jx^ = +Am\ for the FFNS so that both scales be- 
come approximately equal in the limit Q^^mj^; as mentioned the FFNS results are not very sen- 
sitive to this specific choice of the factorization and renormalization scale [48]. Notice that the 
NLO-VFNS predictions for xh (short-dashed curves) are very similar to the ones for {2ej^)^^F2 
(dashed curves) despite the fact that {2e\Y^F2 = (1 + a^Q) ®h+ \asCg ®g, i.e. the ^(a,) quark 
and gluon contributions almost cancel. As expected [48] the discrepancies between the predic- 
tions for xh{x, Q^) in the VFNS and for {2ej^)^^F2{x, Q^) in the FFNS in the relevant kinematic 
region (small x, large Q^) never disappear and can amount to as much as about 30% at very small- 
X, even at = — 1) far above threshold, i.e. 3> ^th- This is due to the fact that here the 
ratio of the threshold energy of the massive subprocess (/*§ — > hh, etc.) and the mass of 
the produced heavy quark V^th/'^/i = 2 is not sufficiently high to exclude significant contribu- 
tions from the threshold region. Even for the lightest heavy quark, h = c, such non-relativistic 
(j8c = \pc\/Ec < 0.9) threshold region contributions to ^2 (-^j Q^) sizeable for < 10^ GeV^ 
due to significant jSc < 0.9 contributions, and the situation becomes worse for h = b [48]. 
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Figure 2.29: The predicted ^-dependencies of I-F2 (x, 2^) in the FFNS and VFNS at some typical fixed values of 
Q^. For the FFNS the renormalization and factorization scales are chosen to be = flp = = Q~ +4m^ with 
mc= 1.3 GeV, and, as usual, = for the VFNS. The NLO-VFNS charm distribution is given by xc{x,Q^) 
as shown by the short-dashed curves. 
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Figure 2.30: As in Fig. 2.28 but for bottom production, i.e.|F^(x, Q^), choosing fl^ = flp = 11^ = +4m| with 
= 4.2 GeV for the FFNS. The short-dashed curves show the NLO-VFNS bottom distribution xb{x, Q^). 
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Figure 2.31: LO predictions for the x-dependencies of the weak charged-current structure function ^F^^{x, Q^) 
for tb production in the FPNS at some typical fixed values of Q^; here the momentum scale is chosen to be 
H^ = IJ,^ = lJ.p = Q^ + {nit + nib)^ with nit = 115 GeV. These predictions are compared with the bottom distribution 
e^ + m?) inthe VFNS where £, =x(l + 



This is in contrast to processes where one of the produced particles is much heavier than the 
other one, like the weak CC contribution W+g^tbio F^'^. Here ^/loi/mb = » 1 and the 
extension of the threshold region where < 0.9, being proportional to mi, / -y/S^ ^ 1 , is strongly 
reduced as compared to ^ = 0.5 in the former case of hh production. Thus the single top 
production rate in W'^ g -^tb is dominated by the (beyond-threshold) relativistic region where 

> 0.9 and therefore is expected to be well approximated by W^b — > t where ^ is an effective 
massless parton within the nucleon. In Fig. 2.30 we compare the LO FFNS^'' predictions for 
jF^^{x,Q^) with the corresponding VFNS ones for S,b{^,Q^ + mj), where the latter refers to 

' 2 

the W^b t transition using the slow rescaling variable [170, 171] ^ = a:(1 + ^) with ntt = 175 
GeV; for F^^{x,Q^) we used ji^ = = fi^ = + {ntt + nth)^. As expected the differences 
between the two schemes are here less pronounced than in the case of cc and bb electroproduction 
in Figs. 2.28 and 2.30. 

'"Notice that the fully massive NLO FFNS QCD corrections to W+ g — » f are unfortunately not available in the literature. 
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2.10 Weak Gauge Boson and Higgs Production at Tevatron and LHC 

The results of Sec. 2.9 indicate that one may resort to the simpler VFNS with its massless h{x, Q^) 
distributions to estimate rather reliably the production rates of heavy quarks, gauge bosons, Higgs 
scalars, etc. at Tevatron and LHC energies, for processes in which the invariant mass of the 
produced system far exceeds the mass of the participating heavy quark flavor. As a next test of 
these VFNS distributions we turn now to hadronic production and present in Fig. 2.31 our 
[4] LO/NLO FFNS/VFNS predictions for o{pp-^W^X) as compared to the data [172-176] and 
to predictions based on the NLO CTEQ6.5 distributions [98]. Although quantitatively slightly 
different, the dominant light quark contributions in the FFNS (wtf— > mJ— > W+, etc.) are due 
to the same subprocesses as in the VFNS, but the relevant heavy quark contributions have been 
calculated via gs{d)^cW^ , gu^bW^ , etc. as compared to cs(<i) — > bu^W^, etc. in the 
VFNS. Here we again expect the VFNS with its effective massless "heavy" quark distributions 
h{x, Q^) to be adequate, since non-relativistic contributions from the threshold region in the FFNS 
are suppressed due to \fffal^cb — ^ 1- Taking into account that the NLO/LO K factor is 
expected [67] to be in the range of 1.2 - 1.3, our LO-FFNS predictions in Fig. 2.31 imply equally 
agreeable NLO predictions as the (massless quark) NLO-VFNS ones shown by the solid and 
dashed-dotted curves in Fig. 2.31. 



■ GJR 
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Figure 2.32: Our (GJR) [4] predictions for the total + production rates at pp colliders. The data have 
been taken from [172-176] and the NLO-VFNS CTEQ6.5 distributions from [98]. The adopted momentum 
scale is /i^ = /i^ = /x| =M^. The scale uncertainty of our NLO GJR predictions, due to varying /i according to 
^Mw < /i < 2Mw, amounts to less than 2% at = 1.96 TeV, for example. 
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total 


107.5 (94.2) 


79.1 (67.9) 


186.5 (162.1) 


~ 111.4(97.0) 


~ 81.2 (69,6) 


~ 192.7 (166.7) 


total 



Table 2.1: NLO(LO) VFNS and FFNS predictions for production at LHC. The uncertainties implied by 
different scale choices are summarized in Eqs. 2.12 and 2.13. The total NLO-FFTSIS rates are obtained by adopting 
an expected [67] AT-factor of 1.2 for the subleading gluon initiated LO rates involving the heavy c and b quarks. 



In Table 2.0 we present our VFNS and FFNS predictions for production at LHC and 

compare the relevant subprocess contributions to G{pp^W'^X) at ^/s = 14 TeV. The light quark 
fusion contributions in the ud and us sector turn out to be somewhat larger in the FFNS than in 
the VFNS which is due to the fact that the u,d,s (and gluon) distributions are evolved for fixed 
nf = 3 in the FFNS. More interesting, however, are the subprocesses involving heavy quarks. 
Here the LO-VFNS is capable to reproduce, to within less than 15%, the LO-FFNS predictions 
based on the gluon-induced fixed-order subprocesses gu^bW, gd^cW and in particular 
on the sizeable CKM non-suppressed gs^cW contribution. As mentioned above, the NLO 
corrections to these latter heavy quark contributions cannot be calculated for the time being. 
However, since these contributions amount to about only 15% of the total FFNS results for 
production (being dominated by the light ud and us fusions in Table 2.0), we can safely employ 
the expected [67] K factor of A'c^ 1.2 for the relevant gd^cW and gs^cW LO contributions 
in Table 2.0 for obtaining the total NLO-FFNS predictions without committing any significant 
error. The resulting total rate for + W~ production at LHC of 192.7±4.7 nb is comparable 
to our NLO-VFNS prediction in Table 2.0 of 186.5 it 4.9 nb where we have added the ±1 (7 
uncertainties implied by our dynamical parton distributions. Using our "standard" FFNS parton 
distributions instead of the dynamical ones for generating the VFNS distributions, the dynamical 
NLO-VFNS prediction slightly increases to 190.7 nb. This latter prediction reduces to 181.0 nb 
when using the smaller scale = \M^. The scale uncertainties of our predictions are defined by 
taking \M\/v < H < 2Mw, using = 80.4 GeV, which gives rise to the upper limits = 2Mw) 
and lower limits (/J = ^Mw) of our predicted cross-sections. In this way we obtain the following 
total uncertainty estimates of our NLO predictions at LHC: 



, r 186.5 ±4.9pdfl^-^|scaienb (VFNS) 
g{pp^W+ + W-+X)= < (2.12) 

I 192.7 ±4.7pdf^3J|3,^,nb (FFNS) 

and, for completeness, at LO 

^ r 162.1 ±3.9pdf;^?i|scaienb (VFNS) 
<7{pp^W+ + W- +X) = i (2.13) 

I 166.7±4.0pdfil^:^|scaienb (FFNS) 
where the subscript "pdf" refers to the lo uncertainties of our parton distribution functions. 
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For comparison, the NLO-VFNS prediction of CTEQ6.5 [98] is 202 nb with an uncertainty of 
8%, taking into account a pdf uncertainty of slightly more than 2<j. Similarly, MRST [84, 85] 
predict about 194 nb. From these results we conclude that, for the time being, the total 
production rate at LHC can be safely predicted within an uncertainty of about 10% irrespective 
of the factorization scheme. 

It is also interesting to study the dependence of the FFNS predictions for the contributions 
to production involving heavy quarks on the chosen scale n as shown in Figs. 2.32 and 
2.34. In these figures we compare the gs^cW initiated production rates in the FFNS with 
the quark fusion cs^W ones in the VFNS and similarly the gd^cW ones with the cd^W 
fusion respectively. These factorization scheme dependencies are rather mild for the LO-FFNS 
predictions, in contrast to the situation for the LO-VFNS predictions which stabiUze, as expected, 
at NLO. The mild n dependence is similar to the situation encountered in tW production [67] via 
the subprocess gb^tW~. 

A similar situation, where the invariant mass of the produced system sizeably exceeds the 
mass of the participating heavy quarks, is encountered in (heavy) Higgs H production accom- 
panied by two heavy fc-quarks, for example. Here H = H^y^; hP, H^,A^ denotes the SM Higgs- 
boson or a Ught scalar hP, a heavy scalar and a pseudoscalar of supersymmetric theories 
with Mh ^100 GeV. In the FFNS the dominant production mechanism starts with the LO sub- 
process gg bbH {qq bbH is marginal), to be compared with the bb fusion subprocess 
in the VFNS starting With bb^ Hat LO. Again, ^/§g/mb = ^"'^J'^^ ~ » 1 in the FFNS 
which indicates that the simpler LO and NLO (NNLO) VFNS bb fusion subprocesses may serve 
as a reasonable effective approximation. Within the scale uncertainties it turns out that the FFNS 
and VFNS predictions at NLO are compatible [177-180]. This result holds for scale choices 
liR^f = {^to\)^/I^ with being considered as a suitable "central" choice in b{x,^p) for 

calculations based on the bb fusion process in the VFNS. It should, however, be mentioned that 
the VFNS tends to somewhat overestimate [177-180] the exact FFNS Higgs-boson production 
rates by about^° 10 to 20%. 

Finally let us note that all our results and comparisons concerning the VFNS hold irrespective 
of the specific parametrizations used for the "heavy" h{x, Q^) distributions. When comparing our 
VFNS distributions [4] with the ones of CTEQ6.5 [98], the relevant ratios vary for lO""^ < x < 0. 1 
at most between 0.9 - 1.1 at LO and NLO; similar results hold for other VFNS distributions, e.g., 
those of [103]. This is illustrated more quantitatively in Fig. 2.34, where we compare our c, b 
and gluon distributions with the ones of CTEQ6 [95] and CTEQ6.5 [98] in the relevant region of 
X < 0.3 at Q^=My^^. The ratios for the light u- and ^/-distributions are even closer to 1 than the 
ones shown in Fig. 2.34, typically between 0.95 and 1.05. 

''Note that the independent variation of iip and /is considered in [177-180] is, as mentioned before, not strictly compatible with 
the utilized parton distributions determined and evolved according to ^.r = Hf- 
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2. Analysis and Results 
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Figure 2.33: The scale (n^ = Hp = jJi^) dependence of the LO-FFNS contribution to the total W+ + pro- 
duction rates due to the subprocesses gs^cW and gd^cW compared to the LO and NLO ones in the VFNS 
due to cs—^W and gd^cW fusion respectively. The results refer to the pp-lMC {■\/s= 14 TeV) and to the 
p^Tevatron (v^= 1.96 TeV) with the latter ones being multiplied by the indicated factor of 10 as indicated. 
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Figure 2.34: As in Fig. 2.32 but for the FFNS subprocess gd—>cW tobe compared with cd - 
The results for the Tevatron {■\/s= 1.96 TeV) are multipUed by a factor of 5 as indicated. 



► Win the VFNS. 
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2.11. Isospin Violations and the NuTeV Anomaly 




Figure 2.35: Comparing our present (GJR-VFNS) dynamical parton distributions generated in the YFNS at 
NLO (MS) with the ones of CTEQ6 [95] and CTEQ6.5 [98] at 2^ = M|,. 

2.11 Isospin Violations and the NuTeV Anomaly 

The NuTeV collaboration has reported [181] a measurement of the Weinberg angle = sin^ Ow 
which is approximately three standard deviations above the world average of other electroweak 
measurements ("anomaly"), i.e. sin^ Ow = 0.2277 ± 0.0013 ± 0.0009 as compared to the world 
average of sin^ Ow = 0.2228 ± 0.0004. Although the NuTeV group has included several uncer- 
tainties in their original analysis, e.g. effects due to a non-isoscalar target, higher-twists, charm 
production, etc.; several sources, spanning from new physics to nuclear effects, have been pro- 
posed to explain this discrepancy (see, for instance, [182-185]). Among them isospin-symmetry 
violating contributions of the parton distributions in the nucleon and effects caused by the strange 
sea asymmetry s ^ s turn out to be relevant. According to the experimental methods used for the 
extraction of the Weinberg angle from measurements of 



,2^V(V)]V/ ^2) 



(2.14) 



cc 

2 



the NuTeV collaboration have presented [182] functionals F[s^,8q;x] (of. Fig. 2.35) which per- 
mit the evaluation of the shifts in the NuTeV result for caused by a momentum distribution 
asymmetry 5q{x,Q^) as: 

a4= / F[s^,dq;x\dq{x,Q^)dx (2.15) 

where the asymmetries are evaluated at (2^ — 10 GeV^ appropriate for the NuTeV experiment. 
We [1, 2] shall approach the issue of isospin violations following Sec. 1.5 within the frame- 
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Figure 2.36: The functionals describing the shift in the NuTeV measurement caused by isospin and strange 
sea asymmetries in the nucleon as presented in [182]. 

work of the dynamical (radiative) parton model and obtain predictions which depend on a single 
free parameter, required by the nonperturbative contributions and chosen to be niq = 10 MeV, 
i.e. of the order of the current quark masses, which are the usual kinematical lower bound for a 
photon emitted by a quark (similar to the electron mass ntg for a photon radiated off an electron 
[186]). We integrate Eq. 1.32 as: 



and similarly for 5m and 5d. To evaluate this equation we use the LO parton distributions (x, q^) 
and q{x,q^) of the dynamical (radiative) parton model from [119]^^. The parton distributions at 
q^ < /^LO' where /x^o ~ ^-^^ GeV^ is the input scale in [1 19], are taken to equal their values at 
the input scale /x^o- resulting valence and sea asymmetries at 2^ = 10 GeV^ are presented in 
Figs. 2.36 and 2.38 respectively. The valence asymmetries are compared with nonperturbative bag 
model results from [185] which arise from an entirely different origin, namely through the mass 
differences dm = md — m„ and dM = M„ — Mp. As can be seen, our radiative QED predictions 
and the bag model estimates are comparable for duy but differ considerably for Sdy. Although 

'^Note that our new dynamical distributions where not available at the time of making our original analysis [2] of the NuTeV 
anomaly. Our results are however very similar to the GRV98 ones (cf Sec. 2.4). 




(2.16) 
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2.11. Isospin Violations and the NuTeV Anomaly 
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Figure 2.37: The isospin violating "majority" Sui, and "minority" 8dy valence quark distributions evaluated at 
Q^ = lOGsV^. Ourradi 
been taken from [185]. 



Q^ = 10 Ge . Our radiative QED predictions are calculated according to Eq. 2 . 1 6 ; the bag model estimates have 
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Figure 2.38: The isospin violating sea distributions 5u and dd at = 10 GeV^ calculated according to Eq. 2.16 
with Uy, dy replaced by u,d respectively. 
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5uv 


5dv 


5u 


8d 


total 


QED 


-0.00071 


-0.00033 


-0.000019 


-0.000023 


-0.0011 


bag 


-0.00065 


-0.00081 






-0.0015 



Table 2.2: The QED corrections to Ay^ evaluated according to Eq. 2.15 using Eq. 2.16. The nonperturbative bag 
model estimates [187] are taken from [185]; different nonperturbative approaches give similar results [185]. 

our method differs somewliat from a similar evaluation [30] of the QED isospin violations, our 
resulting 5qv{x, Q^) turn out to be quite similar as already anticipated there [30]. 

Turning now to the impact of our results on the NuTeV anomaly we present in Table 2. 1 the 
corrections to evaluated according to Eq. 2.15. For comparison, the QED valence isospin 
asymmetries of [30] imply, via Eq. 2.15, contributions to Asy^ similar to ours, namely - 0.00090 
and - 0.00043 due to duy and ddy, respectively. Also shown in Table 2.1 are the additional contri- 
butions to As^f stemming from the nonperturbative hadroiuc bag model calculations. In addition 
we use a nonperturbative estimate [188] for the effects caused by the strange sea asymmetry. Both 
additional contributions are comparable in size to our radiative QED results; the total correction 
becomes: 

2 I 2 I 2 I 2 I 

^■^wltotal — As'^^Iqed + A5(y Ibag + Ai'jy Istrange 

= -0.0011 - 0.0015 - 0.0017 = -0.0043 (2.17) 

Thus the NuTeV measurement of sin^ Ow = 0.2277( 16) would be shifted to sin^ Ow = 0.2234( 16), 
which is in agreement with the standard value 0.2228(4). 
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Summary and Conclusions 



Utilizing recent NC DIS structure function measurements and data on hadronic Drell-Yan dilep- 
ton production and high-px inclusive jet production, we have updated the previous GRV98 [119] 
global fits for the dynamical parton distributions of the nucleon at LO and NLO of perturbative 
QCD and extended our analyses with an appropriate uncertainty analysis for each of our fits. 
The former dynamical distributions are compatible with the new ones at the 1 to 2<j level of the 
newly estimated uncertainties. Furthermore, we have extended our analyses to the NNLO, where 
the high-pT jet data as well as the heavy quark electroproduction data have been disregarded for 
consistency reasons, i.e. since NNLO corrections have not yet been calculated. 

The small-x structure of the dynamical distributions is generated entirely radiatively from 
valence-like, definite positive, input distributions at an optimally chosen input scale /^o < ^ GeV^. 
Our NLO results are stable with respect to the choice of factorization scheme (MS versus DIS) 
and the NNLO predictions are perturbatively stable with respect to the NLO ones and are in agree- 
ment with all present measurements for Q}^1 GeV^. In general, the NNLO corrections imply an 
improved value of x^, typically ;tNNLO — ^-^/^nlo- However, even though the dynamical NNLO 
uncertainties are somewhat smaller than the NLO ones, it turned out that present DIS precision 
data are still not sufficiently accurate to distinguish between NLO results and the minute NNLO 
effects of a few percent. The stability of these dynamical results guarantees reliable calculations 
of cross-sections for, e.g., heavy quark, W^, Z^, and high-px jet production at hadron colliders 
like Tevatron and in particular LHC. 

Our dynamical distributions have been compared with conventional "standard" ones obtained 
from non-valence-like positive gluon and sea input distributions at some arbitrarily chosen higher 
input scale > 1 GeV^. For this purpose we have also performed common "standard" fits in ex- 
actly the same conditions as the dynamical ones but assuming I^q=2 GeV^. Notice that, contrary 
to the dynamical approach, the finite small-.^: behavior of the input gluon and sea distributions is 
here fitted, and not radiatively generated by QCD evolutions. For this reason as well as because 
of the sizably different evolution distances, i.e. sizably different Hq, the uncertainties of these 
less constrained "standard" distributions are larger than those of their dynamical counterparts, in 
particular in the small-x region. 

The more constrained ansatz of the dynamical model implies a certain predictive power, in 
particular in the present experimentally unexplored extremely small-x region relevant for eval- 
uating ultrahigh energy neutrino-nucleon cross-sections in astrophysical applications. Here we 
provide predictions down to xc± I0~^ at the weak scale Q^ = My^^ as required [160-162] for the 
highest cosmic neutrino energies of 10^^ GeV. These predictions are strongly constrained within 
the dynamical parton model and are entirely of QCD-dynamical origin in the very small-A: region. 
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Since previous predictions [1 16, 1 17, 122] for the small-x region based on the dynamical parton 
model and the data available at the time were subsequently confirmed [120, 121] at HERA, the 
presently available very precise small-x data [106, 123-126] utilized here allows us to be quite 
confident about the reliability of our improved small-x predictions within the framework of the 
successful dynamical parton model. 

Our predictions for the longitudinal structure and results for the "reduced" DIS cross- 
section Gr are in agreement with all HERA data, including the most recent HI measurements, 
in particular in the small-x region and down to = 1 GeV^. The dynamical NLO/NNLO 
predictions for Fl are positive throughout the whole kinematic region considered and become 
perturbatively stable already at (2^ = 2 — 3 GeV^, where future precision measurements could 
even delineate NNLO effects in the very small-x region, where the NNLO uncertainty bands are 
smaller than the NLO ones. This is in contrast to the common "standard" approach, however 
instabilities and differences are here less distinguishable due to the much larger uncertainties. 

In order to analyze the role of heavy quark flavors in high-energy colliders, we have generated 
radiatively sets of (zero-mass) VFNS parton distributions based on our FFNS LO and NLO dy- 
namical parton distributions. We have confronted these VFNS and FFNS predictions in situations 
where the invariant mass of the produced system {hh, tb, cW, tW, Higgs-bosons, etc.) does not 
exceed or exceeds by far the mass of the participating heavy flavor. In the former case the VFNS 
predictions deviate from the FFNS ones by up to about 30% even at ml while in the latter 
case these deviations are appreciably smaller, within about 10%, which is within the margins of 
renormaUzation and factorization scale uncertainties. Our results indicate that the simpler VFNS 
with its effective treatment of heavy quarks (c, b) as massless partons can be employed for calcu- 
lating processes where the invariant mass of the produced system is sizeably larger than the mass 
of the participating heavy quark flavor. Taking into account the uncertainties of parton distribu- 
tions and those due to scale choices, the total production rate at LHC can be predicted within 
an uncertainty of about 10%, which becomes significantly smaller at the Tevatron. Similarly the 
Higgs production rates at LHC are predicted with an uncertainty of 10 - 20% where the VFNS 
production rates exceed the FFNS ones by about 20% at LHC. 

Finally, we have evaluated the modifications to the standard isospin symmetric parton distri- 
butions due to QED ^(a) photon bremsstrahlung corrections within the dynamical model and 
compared them with nonperturbative bag model calculations, where the violation of isospin sym- 
metry arises from entirely different (hadronic) sources, predominantly through quark and target 
mass differences. Taken together these two isospin violating effects and effects caused by a 
strangeness asymmetry (s^s) in the nucleon the discrepancy between the large result for sin^ 9w 
as derived from deep inelastic v(v)A'^ data ("NuTeV anomaly") and the world average of other 
measurements is entirely removed. 

FORTRAN codes (grids) containing our dynamical parton distribution functions along with 
the uncertainties of the most relevant sets can be obtained on request^^ or directly downloaded 
from http://dooin.physik.uni-dortinund.de/pdfserver; they are also included in the Les 
Houches Accord PDF Interface [189]. We hope that they will be of utility for the analysis and 
understanding of the (new) results which will be coming soon in High-Energy/Particle Physics. 



'pjimenez@het.physik.um-dortmund.de 
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A. Tables of Numerical Values 



We present here the numerical values obtained in our different fits and error analyses. Tab. A.l 
represent a summary of our results; there we include the central values obtained for the parameters 
in our fits, as{M^) and the parameters of the input distributions, and their errors. It also contains 
the minimum values of obtained for the complete fits as well as for different combinations 
of data which we find interesting to quote separately. We did not find that useful to include 
the obtained for all individual data sets; nevertheless the interested reader can obtain them 
on request^ In Tab. A.l we quote the values obtained for the overall normalization factors for 
each individual data set in our different fits, as already explained (cf. Sec. 2.2), this factors were 
allowed to float within the experimental normalization uncertainty of each data set. 

Finally, in Tables A.3 to A.9 we give all the information necessary for uncertainty studies (cf. 
Sec. 2.3). Each of these tables contains the eigenvalues of the hessian matrix and the eigenvector 
basis sets for one of our fits; the central values are included again for comparison. Recall that 
with these values it is possible to calculate the uncertainty estimated for any quantity as well as its 
derivatives at the central values of the parameters in the linear approximation. The reader inter- 
ested may also reconstruct the coefficients (/?/), which give information on the correlations of 
the different parameters, or even the hessian matrix and therefore the (approximated) x^ profiles 
of our fits. 
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Table A.2: Overall normalization factors (in percent) multiplying the data as obtained in our different fits. Obvi- 
ously, a quoted value of Aget means that each data point of the set is to be multiplied by a factor A^set = 1 + O.OlAset 
in order to obtain the 'x} values quoted in Tab. A.l. 
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B. Splitting Functions 



We present here the expressions for the sphtting functions which have been implemented in our 
code for the solutions of the RGE (cf. Sec. 1.4). Together with the well-known one-loop splitting 
functions [191, 192], and the two-loop ones [37, 190, 193, 194], the expressions for the 3-loop 
splitting functions [195, 196] are necessary for the evolution of parton densities up to NNLO. 
For these latter quantities we use the ;ic-space parametrizations presented also in [195, 196], and 
transform them to w-space using, e.g., [197]. For the rip contribution in the non-singlet sector 
we use, however, the exact «-space expression as presented in [198]. 

These quantities can be expressed in terms of simple functions and harmonic sums [190, 197]; 
the functions that we use are defined in terms of integer n values as [190]: 

'■=^ _ . (B.l) 

5Ui)=5.(int(i))=2'«-i£l±l^. 
^ ^ (=1 ' 

Their analytic continuations to complex-n space can be expressed in terms of logarithmic deriva- 
tives of Euler's F-function. For the first harmonic sums they are given by: 

Si{n) = rE + V^^\n+l), 7^ = 0.577216... , 

52(n) = C(2)-V^W(n+l), C(2) = ^, (B.2) 

53(n) = C(3)-^V^(2H«+l)> C(3) = 1-202057... , 

where \l/^'\n) = '^''^Ji+i"^ ■ Since in principle S'^{^) and S{n) are defined for even (odd) n, there 
are two analytic continuations, which we denote + (-); they are given by: 

S.W = -if (3)± (aw _ ia (^(^) _ ^(D) y?M) . 

For the calculation of the above integral we use in practice the following parametrization [1 16]: 



i 



^ n-lLilW ^ ^ 1.01 _ 0.846 _L 1-155 _ 1.074 _l 0.55 

' 1 +x 



^yyW-i ^'^V'V I l.Ul U.846 I 1.155 1.U74 , U.55 (R A\ 
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Finally, we list the expressions used for the splitting functions in the remaining of the appendix : 



p(0) 



n(n+l) 



(B.5) 



(0) _ 8 / n^+n+2 \ 



,(0) 



3 |(n-l)«(«+l) J' -^gg 



^{T + n(h) + (n+lfci) -45i(n)} - In, (B.6) 



4^'^ = -^js^iW;;^ + 8 [25i(n) - ^] [52(n) + 1252(n) 

+325±(n)-45^^i(f)-|-4 



-4<! +^5i(n) 



25i(n) 



. «3(„+l)3 
1 



t8 



n{n+l) 



252(«)-5^(f 



44 

3 



-165±(n) + 25^^^(^)-^-^ 
-fnf|-f5i(n) + ^52(n) + f + | 



15 ln*+236n3 +88n^+3n+ 1 8 



n3(n+l)^ 
ll«^+5n-3 



+ 4 



52(«) 

2n^+2«+l 



(B.7) 



p(i) 

^PS 



16 ) irr' +32n*+A9r? +38n^+28>i+8 
3 '"l (n-l)n3(„+l)3(„+2)2 



(B.8) 



p(i) 
^qg 



+6?if |2 



S2{n)-S'^^^{l)-S\{n) 



rp-+n+2 



+2 



n(n+l)(n+2) 

>i^+6«^+15n'^+25?i^+36n^+85«4+l28«3+l04n^+64n+16 



+ 85i(n) 



2«+3 
(n+l)^(n+2)^ 



(«-l)n3(n+l)^(n+2)^ 



[25?(«)-252(n) + 5] 

lln^+26»3+15n^+8ra+4 1 
n3(«+l)3(„+2) J 



«^+n+2 



n(n+l)(n+2) 



-4 



5i(«) 



(B.9) 



^g? = + f [-25?(«) + 105i(n) -252(n)] 



rt+n+2 
{n-\)n{n+\) 
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12>i^+30«^ +43n'* +28^3 -w^ - 1 2W-4 



+ 16 



(n-l)n3(n+l)3 

52(n) + 52(n)-5' 



n^+n+2 
{n-V)n{n+\) 



-5i(n) 



17«'*+41«^-22«-12 



3(n-l)V(n+l) 

10 V +621n^ + 1400«'' + 1 678n^ +695n5 - 1 03 1 n"* - 1 304n3 _ 1 52«2+432n+ 144 
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'Note that the numerical values of the QCD color factors are explicitly given, in particular for the 3-loop functions they are 
included in the parametrizations. 
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